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Abstract 


\ 

'"'An  imaging  system's  object  field  to  image  field  trans¬ 
formation  operation  is  usually  described  in  the  (spatial) 
frequency  domain  using  the  system's  transfer  function.  The 
limitation  of  this  description  is  that  the  imaging  system 
must  be  space-invariant.  In  contrast,  an  object  field  to 
image  field  transformation  can  be  described  in  the  spatial 
domain  without  requiring  that  the  system  be  space-invariant. 

For  the  frequency  domain,  a  summary  of  free-space  and 
turbulent  transfer  functions  is  presented.  For  the  spatial 
domain,  the  Normal  Mode  Approach  to  imaging  is  described 
followed  by  a  summary  of  spatial  eigenvalue  distributions 
and  degrees  of  freedom  expressions.  The  effect  of  additive 
background  noise  on  the  useable  degrees  of  freedom  of  an 
image  field  is  studied.  The  spatial  domain  and  frequency 
domain  are  shown  to  be  related  when  certain  conditions  are 
satisfied,  one  of  which  is  space-invariance. 

The  operation  of  two  ideal,  adaptive  imaging  receivers 
(the  Channel-Matched  Filter  and  Multiplicative-Phase  Receivers 
is  described  and  their  imaging  performance  is  compared  using 
turbulent  coherent  transfer  functions  and  minimum  average 
integrated  mean  square  error  expressions. Av 
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I  Introduction 


Background 

Imaging  systems  are  typically  composed  of  a  channel 
through  which  an  object  field  is  propagated  and  a  receiver  at 
the  channel  output  for  imaging  the  object  field.  The  manner 
in  which  these  systems  transform  an  object  field  into  an  im¬ 
age  field  is  'commonly  described  in  the  (spatial)  frequency 
domain  by  using  transfer  functions.  A  transfer  function  indi¬ 
cates  how  the  Fourier  (sinusoidal)  components  of  an  object 
field  are  transformed  into  the  corresponding  components  of 
an  image  field.  If  an  object  field  is  spatially  coherent  the 
function  is  called  a  Coherent  Transfer  Function  (CTF).  If  an 
object  field  is  spatially  incoherent  the  function  is  called 
an  Optical  Transfer  Function  (OTF )  and  its  modulus  is  called 
the  Modulation  Transfer  Function  (Ref  7:114). 

In  the  frequency  domain,  the  resolution  capability  of  an 
imaging  system  is  sometimes  defined  in  terms  of  the  highest 
nonzero  spatial  frequency  (spatial  frequency  cutoff)  of  a 
CTF  or  MTF.  The  higher  the  spatial  frequency  cutoff  the 
greater  the  resolution  capability  of  an  imaging  system.  For 
propagation  of  an  object  field  through  free-space  the  spatial 
frequency  cutoff  is  due  to  the  diffraction  effects  of  the 
finite  apertures  of  an  imaging  receiver.  The  free-space  cut¬ 
off  is  the  highest  resolution  attainable  by  an  imaging  sys¬ 
tem. 

When  the  lack  of  space-invariance  in  an  imaging  system 
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prevents  the  use  of  transfer  functions,  it  is  useful  to  be 
able  to  describe  imaging  in  another  domain  called  the  spatial 
domain.  In  the  spatial  domain,  an  object  field  to  image  field 
transformation  is  described  in  terms  of  object  field  and  image 
field  spatial  modes  and  spatial  eigenvalues.  A  spatial  eigen¬ 
value  indicates  the  power  loss  associated  with  a  single  object 
field  spatial  mode  to  image  field  spatial  mode  transformation. 
When  the  spatial  eigenvalues  are  ordered  in  magnitude  (highest 
to  lowest),  they  form  a  distribution  that  characterizes  the 
modal  transformation  of  an  imaging  system. 

In  the  spatial  domain,  the  resolution  capability  of  an 
imaging  system  is  usually  defined  as  the  number  of  nonzero 
eigenvalues  in  the  spatial  eigenvalue  distribution.  This  num¬ 
ber  is  called  the  degrees  of  freedom  (DOF)  of  an  image  field. 
The  greater  the  number  of  DOF  of  an  image  field,  the  greater 
is  the  resolution  capability  of  an  imaging  system.  The  DOF 
of  an  image  field  depend  on  the  geometry  of  an  imaging  system, 
the  wavelength  of  the  propagated  object  field  and  the  propaga¬ 
tion  medium. 

Another  spatial  domain  resolution  measure  is  the  inte¬ 
grated  mean  square  error  (IMSE)  of  the  image  field.  The  error 
is  the  difference  between  the  image  field  produced  for  free- 
space  propagation  and  either  the  image  field  produced  for 
propagation  in  other  media  or  for  additive  background  noise. 
The  IMSE  of  the  image  field  is  expressible  in  terms  of  the 
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DOF  of  the  image  field 


Objectives 

This  thesis  has  three  objectives.  The  first  objective  is 
to  describe  the  Normal  Mode  Approach  to  spatial  domain  imag¬ 
ing  and  to  present  spatial  eigenvalue  distributions  and  DOF 
expressions  associated  with  that  approach.  The  second  objec¬ 
tive  is  to  show  how,  and  under  what  conditions,  the  spatial 
domain  is  related  to  the  frequency  domain.  The  third  objec¬ 
tive  is  to  compare  two  ideal  adaptive  imaging  receivers  using 
their  turbulent  CTFs  and  the  minimum  average  IMSEs  of  the  image 
fields  they  generate.  The  receivers,  to  be  described  later, 
are  the  Channel-Matched  Filter  (CMF)  Receiver  and  the 
Multiplicative-Phase  (MP)  Receiver. 

Organization 

Chapter  II  is  a  summary  of  the  CTFs  and  MTFs  for  an  ob¬ 
ject  field  propagated  through  either  free-space  or  the  turbu¬ 
lent  atmosphere. 

Chapter  III  is  a  description  of  the  Normal  Mode  Approach 
to  spatial  domain  imaging,  a  summary  of  spatial  eigenvalue 
distributions  and  DOF  expressions,  and  an  examination  of  the 
effect  additive  background  noise  has  on  the  useable  DOF  of  an 
image  field.  In  addition,  Chapter  III  shows  how  the  spatial 
domain  and  the  frequency  domain  are  related. 

Chapter  IV  is  a  description  of  the  CMF  and  MP  Receivers 


and  a  derivation  of  turbulent  CTFs  and  image  field,  minimum 
average  IMSE  expressions  for  both  receivers.  This  chapter 
is  also  a  comparison  of  the  resolution  capabilities  of  the 

receivers. 

Chapter  V  is  a  summary  of  conclusions  and  a  discussion 
of  areas  requiring  further  study . 


II  Frequency  Domain  Imaging-Transfer  Functions 


The  purpose  of  this  chapter  is  to  present,  as  background 
material,  the  coherent  transfer  functions  (CTFs)  and  modula¬ 
tion  transfer  functions  (MTFs)  of  an  object  field  propagated 
through  either  free-space  or  the  turbulent  atmosphere.  The 
first  section  of  this  chapter  defines  the  imaging  system 
geometry  used-  in  calculating  the  CTFs  and  MTFs.  The  next 
section  is  a  summary  of  free-space  CTFs  and  MTFs,  and  corre¬ 
sponding  spatial  frequency  cutoffs,  for  the  defined  geometry. 
The  final  section  is  a  summary  of  turbulent  CTFs  and  MTFs. 

General 

Figure  1  shows  an  imaging  system  composed  of  an  object 
field  plane  P.^,  receiving  plane  P2,and  an  image  field  plane 
P3.  P2  has  a  circular  aperture  R2  of  diameter  d2.  The  region 
between  P.^  and  P2  is  called  the  propagation  channel.  The 
channel  may  be  either  free-space  or  the  turbulent  atmosphere. 
The  region  between  P2  and  P3  is  the  imaging  receiver.  This 
region  could  include  optical  elements  but  for  simplicity  it  is 
assumed  to  be  a  free-space  channel  without  such  elements. 

The  receiving  plane  P2  is  therefore  both  the  input  to  the 
receiver  and  the  output  of  the  propagation  channel.  The  imag¬ 
ing  system  of  Figure  1  is  assumed  to  be  linear  and  space  in¬ 
variant.  This  assumption  must  be  made  if  transfer  functions 
are  to  be  used  to  characterize  the  object  field  to  image  field 
transformation  of  the  imaging  system  (Ref  7:19). 
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Figure  1  -  Geometry  for  Calculation  of  Transfer  Functions 

Free-Space 

Assume  that  the  propagation  channel  of  Figure  1  is  free- 
space.  Associated  with  the  circular  receiving  aperture  R2 
of  Figure  1  is  a  free-space  pupil  function  called  P(?2): 

P(?2)  =  circ(2j?2 |/d2)  (1) 

where  , 

1  ;  x  ^ 1/2 

circ(x)  =  < 


0 


otherwise 


(2) 


Goodman  has  shown  (Ref  7:111)  that  for  a  coherent  object 


field  in  P^,  the  free-space  CTF,  called  ^q-COH^^’  is 

to-COH(?)  "  P(>'dil?l)  <3> 

where  X  is  the  optical  wavelength  of  the  object  field,  d^  is 
the  length  of  the  propagation  channel,  and  j?|  is  a  spatial 
frequency  magnitude.  Combining  Eqs  (1)  and  (3)  yields 

Zo_coh(^  =  circ(2  Adjf  |/d2)  (4) 


From  Eq  (4),  the  CTF  spatial  frequency  cutoff,  called 


?o|  *  d2/2>1d1 


Goodman  has  shown  (Ref  7:116)  that  for  an  incoherent  ob¬ 
ject  field  in  P^,  the  free-space  MTF,  called  7Tq(?)»  is  the 
convolution  of  the  pupil  function  of  Eq  (1)  with  itself. 
Completing  this  convolution  yields 


cos  _1(Adi|f  |/d2)  -  (Ad±|?|/d2) 


T0(?)  *  x  -yfi  -  (Adjf|/d2)2  1/7T  ;  Xdi|?|^d2 


0  ;  otherwise 


(6) 


t 


t 


From  Eq  (6),  the  MTF  spatial  frequency  cutoff,  called 
is 


d2/Xd. 


(7) 


A  comparison  of  Eqs  (5)  and  (7)  indicates  that  the  spa¬ 
tial  frequency  cutoff  for  an  incoherent  object  field  is  twice 
the  cutoff  for  a  coherent  object  field.  This  does  not  imply, 
however,  that  incoherent  object  fields  yield  "better"  image 
fields  than  coherent  object  fields.  As  Goodman  points  out 
(Ref  7:125),  the  two  cutoff  frequencies  describe  resolution 
for  different  image  quantities  and  are  therefore  not  compa¬ 
rable.  Equation  (5)  is  the  spatial  frequency  cutoff  for  a 
function  which  produces  an  image  field  amplitude  distribution 
while  Eq  (7)  is  the  cutoff  for  a  function  which  produces  an 
image  field  intensity  distribution  (power  density). 

Turbulent  Atmosphere 

Assume  now  that  the  propagation  channel  of  Figure  1  is 
the  turbulent  atmosphere.  An  object  field  propagating  through 
the  turbulent  atmosphere  is  effected  by  multiplicative  log- 
amplitude,  X(r2»^i^  and  phase,  0(?2,r1),  fluctuations. 
X(r2,r1)  and  (f){r^,r^)  are  Gaussian  random  variables 
(Ref  17:209).  Their  first  and  second  order  statistics  are 
given  in  Appendix  A.  Since  X(r2,r^)  and  0(r2,?1)  are  in 
general  functions  of  both  the  object  field  plane  and  the  re¬ 
ceiving  plane,  they  violate  the  space  invariance  required 


for  obtaining  transfer  functions.  There  is  a  condition,  how¬ 
ever,  under  which  X(r2,r^)  and  Q(r2,r^)  can  be  expressed  as 
functions  of  only  the  receiving  plane  so  as  not  to  violate 
the  space  invariance  of  the  imaging  system.  The  condition  is 
that  the  maximum  spatial  extent  of  an  object  field  in  P 1  of 
Figure  1  be  less  than  an  atmospheric  coherence  length  rQ 
(Ref  15:462).  tor  a  spherical-wave  object  field,  rQ  is  de¬ 
fined  as  (Ref  8:550) 


r 


o 


0.423 


(8) 


where  k  =  2  717  X  is  the  wave  number  of  the  object  field,  dQ 

2 

is  the  length  of  the  propagation  channel,  and  (z)  is  the 
refractive-index  structure  parameter  of  the  turbulent  atmos¬ 
phere  (Ref  9:1526).  For  X  «  0.55,4  and  a  vertical  propaga¬ 
tion  path,  rQ  is  approximately  11  cm  (Ref  5:2622).  When  the 
object  field  meets  the  above  condition  it  is  said  to  lie  with¬ 
in  an  isoplanatic  patch.  Since  most  object  fields  are  spa¬ 
tially  larger  than  an  isoplanatic  patch,  isoplanatism  is  a 
severe  (but  necessary)  restriction  for  obtaining  a  turbulent 
CTF. 

When  isoplanatism  is  a  valid  assumption,  there  is  asso¬ 
ciated  with  the  circular  receiving  aperture  R2  of  Figure  1 
a  "turbulent"  pupil  function  called  P(?2): 

P(r2)  *  P ( rg ) exp  [X(  |  ?2  1 )  +  j  |  ?2 1  )]  (9) 
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where  P(r2)  is  defined  by  Eq  (1).  Goodman  has  shown 

(Ref  7:121)  that  for  a  coherent  object  field  in  of  Figure  1 

the  turbulent  CTF,  called  is 

Zcm(f)  =  P(Ad.|f|)  (10) 

Combining  Eqs  (3),  (9)  and  (10)  yields 

^COH^)  =  ^0-C0H^^exp  ^  di  I  +  jflcAdJ?!)] 

(11) 

The  ensemble  average  of  Eq  (11)  with  respect  to  the 
random  variables  X(Adi|^|)  and  (^(Ad.jjf])  is 

<^C0H^*'>  =  ^O-COH^  <exP  [X(  A  di  jf|  ) 

+  j  0(  AdJ?|)]  >  (12) 

It  is  shown  in  Appendix  A  that  Eq  (12)  is  equivalent  to 

<^'C0H^  >  =  ^O-COH^^ ex^  [~ ^2] 

2  2  2  2  2 
where  V"  =  ^  ,  and  and  are  the  variances 

of  X(r^)  and  0(?2)  respectively.  Following  Fried's  example 

(Ref  4:1375),  the  exponential  in  Eq  (13)  can  be  thought  of  as 

an  atmospheric  CTF,  called  ^a-COH^*  Equation  (13)  is  then 
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(14) 


<^COH(?)>  =  ^0-C0H(^  Ta-COH(?) 

For  a  spherical  wave,  Greenwood  has  indicated  (Ref  8:550) 
that  \r2  of  Eq  (13)  is 

5 

ST2  »  4.93  (Lo/rQ)3  (15) 

where  rQ  is  the  atmospheric  coherence  length  defined  by  Eq  (8) 
and  Lq  is  the  outer  scale  of  the  turbulent  atmosphere.  A 
typical  value  for  Lq  is  100  meters  or  one-fifth  the  height 
above  the  ground,  whichever  is  less  (Ref  9:1527).  Substitut¬ 
ing  Eq  (15)  into  Eq  (13)  gives 

<'^'C0H^'>  =  ^0-C0H^f^exP[“2*47^Lo/,ro^  ] 

Since  the  exponential  in  Eq  (16)  is  not  spatial  frequency 

dependent,  the  spatial  frequency  cutoff  of  the  turbulent 

CTF  is  the  same  as  the  spatial  frequency  cutoff  of  the  free- 

space  CTF.  However,  the  amplitude  of  the  turbulent  CTF  is, 

for  all  spatial  frequencies,  smaller  than  the  amplitude  of 

the  free-space  CTF  since  L-  is  larger  than  rQ.  The  amplitude 

attenuation  is  usually  due  to  turbulence-induced  phase  fluc- 

2 

tuations  rather  than  to  log-amplitude  fluctuations  since  ^ 

p 

is  usually  much  larger  than  ^  (Ref  8:550). 

In  reference  4,  Fried  has  derived  expressions  for  long 
and  short-exposure  turbulent  MTFs.  Long-exposure  means  that 
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the  object  field  propagated  through  the  turbulent  atmosphere 
is  viewed  over  a  long  enough  time  such  that  the  turbulent  MTF 
is  equal  to  its  ensemble  average.  Short-exposure  means  that 
the  propagated  object  field  is  not  viewed  long  enough  to 
justify  equating  the  turbulent  MTF  to  its  ensemble  average. 

For  short  exposures,  Fried  has  considered  both  near-field  and 
far-field  object  fields.  V/ith  reference  to  Figure  1,  near¬ 
field  means  d2  »  while  far-field  means  a2  <<  Vd  X  . 

Only  the  near-field  case  is  presented  here. 

Fried's  ensemble-average,  long-exposure,  turbulent  MTF, 
called  <rLE(f»,  is 

<tL£(f)>  =  Xq  (  f )  exp  ^3 . 44  (  X  d^  J  f  J  /rQ ) 3  ]  (17) 

vrhere  d„.  is  defined  in  Figure  1  and  rQ  is  defined  by  Eq  (8). 
The  average  of  Eq  (17)  is  an  average  over  many  different  expo¬ 
sures.  Although  the  exponential  in  Eq  (17)  is  similar  in 
form  to  the  exponential  in  Eq  (16),  their  effects  are  differ¬ 
ent.  In  Eq  (16)  the  exponential  represents  a  constant  ampli¬ 
tude  attenuation  for  all  spatial  frequencies  while  in  Eq  (17) 
the  exponential  is  a  spatial  frequency  dependent  attenuation. 

Fried's  ensemble-average,  near-field,  short-exposure, 
turbulent  MTF,  called  <Kp  TSE(?)> ,  is 

^NF^SE^^  =  V?)exp 


5 

<  3.44(Xdi|?|/ro)3 
► 
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X 


(18) 


’ 

[l  -  (\dj?|/d2)2  ]  > 

i 

v/here  is  the  diameter  of  the  circular  receiving  aperture 
of  Figure  1.  An  examination  of  Eq  (18)  indicates  that 
<T>JF  TSE(^)>  approaches  the  free-space  MTF,  as  1^1 

approaches  d2/Xd..  From  Eq  (7),  d2  /X  d^^  is  the  spatial  fre¬ 
quency  cutoff  of  the  free-space  MTF. 

All  of  the  CTF  and  MTF  expressions  in  this  chapter  are 
based  on  the  assumption  that  the  imaging  system  of  Figure  1 
is  linear  and  space-invariant.  When  space- invariance  is  not 
a  valid  assumption,  CTF  and  MTF  expressions  are  meaningless. 
The  loss  of  these  expressions  makes  a  frequency  domain  de¬ 
scription  of  an  object  field  to  image  field  transformation 
very  difficult.  Fortunately  this  transformation  can  be  de¬ 
scribed  in  a  different  domain  without  having  to  assume  space- 
invariance.  The  other  domain  is  called  the  spatial  domain 
and  it  is  the  subject  of  the  next  chapter. 


13 
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Ill  Spatial  Domain  Imaging 


This  chapter  has  four  purposes.  The  first  purpose  is 
to  explain  the  Normal  Mode  Approach  to  spatial  domain  imag¬ 
ing.  The  second  is  to  present  spatial  eigenvalue  distribu¬ 
tions  for  propagation  of  coherent  or  incoherent  object  fields 
through  either  free-space  or  the  turbulent  atmosphere  and  to 
give  expressions  for  the  degrees  of  freedom  (DOF)  of  the  re¬ 
sulting  image  fields.  The  third  purpose  is  to  determine  the 
effect  additive  background  noise  has  on  the  DOF  of  an  image 
field.  The  final  purpose  is  to  show  how,  and  under  what  con¬ 
ditions,  the  spatial  domain  is  related  to  the  frequency  domain. 
The  sections  of  this  chapter  are  arranged  in  the  same  order 
as  the  purposes  above. 

Normal  Mode  Approach 

Figure  2  shows  the  imaging  system  geometry  to  be  used  in 
explaining  the  Normal  Mode  Approach  to  spatial  domain  imaging. 
P.^  is  an  object  field  plane  and  P2  is  an  image  field  plane. 

The  region  between  P^  and  P2  is  called  the  propagation  chan¬ 
nel.  Therefore,  P2  is  not  only  an  image  field  plane  but  also 
the  output  plane  of  the  propagation  channel.  In  later  sec¬ 
tions  of  this  chapter  the  channel  is  either  free-space  or 
the  turbulent  atmosphere,  but  for  now  it  is  a  general  channel 
with  a  general  impulse  response  called  h21^2'^l^*  The  pur¬ 
pose  of  aperture  of  Figure  2  is  to  define  the  spatial  ex¬ 
tent  of  the  object  field  as  it  enters  the  propagation  channel. 


14 


Aperture  R2  serves  the  same  purpose  for  the  image  field  at 
the  output  of  the  propagation  channel.  In  this  section,  the 
shapes  of  and  R2  are  arbitrary  but  in  later  sections  they 
are  slit  or  circular  apertures.  Since  in  the  Normal  Mode 
Approach  the  spatial  extents  of  object  and  image  fields  are 
defined  in  their  respective  planes,  a  receiving  plane  and  a 
receiving  aperture  are  not  required  in  Figure  2  as  they  are 
in  Figure  1.  The  imaging  system  of  Figure  2  is  linear  but 
not  necessarily  space-invariant. 


Figure  2  -  Geometry  for  Normal  Mode  Imaging 


In  the  Normal  Mode  Approach  an  object  field  EjCr^  in  R^ 
of  Figure  2  is  expressed  as  a  weighted  sum  of  a  complete  ortho- 
normal  (CON)  set  of  spatial  eigenfunctions  (spatial  modes)  as 
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shovTi  below: 


oo 

EI(?l)  =  Yja i  (19) 

i=l 

where  the  0^(r^)  are  a  CON  set  of  spatial  eigenfunctions 
defined  in  R^  and  the  .  object  field  weighting  coeffi¬ 

cients  defined  by 


J 

R, 


El(?l) 


<"i 


(r1)d?1 


(20) 


Likewise,  an  image  field  EQ(r2)  in  of  Figure  2  is  ex¬ 
pressed  as  a  weighted  sum  of  a  COM  set  of  spatial  eigenfunc¬ 
tions  (spatial  modes)  as  shown  below: 


oo 

EQ(r2)  =  (21) 

i  =1 

where  the  X^(r2)  are  a  CON  set  of  spatial  eigenfunctions  de¬ 
fined  in  R2  and  the  bi  are  image  field  weighting  coefficients 
defined  by 


(r2)dr 


2 


(22) 


In  addition  to  Eq  (21),  the  image  field  is  expressed  in 
terms  of  the  object  field  and  the  general  impulse  response  of 
the  propagation  channel  by  using  the  Huygens-Fresnel  Principle 
(Ref  2:Chap  6): 
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(23) 


Eo(?2) 


■  /EI(?1 


)h2i(r2,ri)dri 


The  Normal  Mode  Approach  describes  an  object  field  to 
image  field  transformation  by  describing  how  much  of  the 
power  of  each  of  the  spatial  modes  of  the  object  field  is 
transmitted  through  the  propagation  channel  to  the  corre¬ 
sponding  modes  of  the  image  field.  On  a  per  mode  basis  the 
fraction  of  power  transmitted  is  given  by  the  ratio  of 
jbjJ2  to  Ja^l^.  This  ratio  is  defined  as 

MVM8  -  K  <24> 

where  is  called  the  ith  power  spatial  eigenvalue  of  the 
propagation  channel.  In  the  remainder  of  this  thesis  X.^  is 
referred  to  as  simply  the  ith  spatial  eigenvalue.  The  square 
root  of  the  ith  spatial  eigenvalue  is  related  to  the  ith  spa¬ 
tial  modes  of  the  object  and  image  fields,  and  to  the  general 
impulse  response  of  the  propagation  channel,  by 


V\  ^(rg)  =  ^  ^i(?l)h21(?2,?l)d?l 


Equation  (25)  is  derived  in  Appendix  B.  In  Appendix  B 
it  is  shown  that  one  way  Eq  (25)  is  satisfied  is  for 

\  01(?')  =  /k(?1*,?1)  01(P1)dr1  ;  r1'€R1  (26) 
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where 


K(?l ' ,?!>  =  J  h21(r2,T1)h2l*{r2,r1')dr2  ; 

2  (27) 

In  Eqs  (26)  and  (27),  K(r^' ,r^)  is  called  an  imaging 
kernel.  Equation  (26)  is  called  a  Fredholm  Equation.  For 
free-space  propagation,  Slepian  and  Pollack  have  shown 
(Ref  16:57)  that  the  spatial  eigenfunctions  0  i(r1)  that 
satisfy  Eq  (26)  are  called  prolate  spheroidal  wave functions. 
Using  these  wavefunctions,  Eq  (26)  can  be  solved  for  free- 
space  spatial  eigenvalues.  The  eigenvalues  can  be  ordered 
such  that  X,  ^  •  The  ordered  spatial  eigenvalues 

form  a  distribution  that  charset  jrizes  the  modal  transforma¬ 
tion  of  an  imaging  system.  The  number  of  nonzero  spatial 
eigenvalues  in  the  spatial  eigenvalue  distribution  is  called 
the  degrees  of  freedom  (DOF)  of  an  image  field.  The  greater 
the  number  of  DOF  of  an  image  field,  the  greater  is  the  reso¬ 
lution  capability  of  sin  imaging  system. 

For  propagation  through  the  turbulent  atmosphere,  the 
imaging  kernel  and  spatial  eigenfunctions  are  random  varia¬ 
bles,  Since  the  spatial  eigenfunctions  $^(r^)  that  satisfy 
Eq  (26)  are  not  known,  Eq  (26)  cannot  be  solved  for  turbulence 
spatial  eigenvalues.  Therefore,  turbulence  spatial  eigen¬ 
values  are  also  random  variables.  Shapiro  has  shown 
(Ref  13:2616),  however,  that  turbulence  spatial  eigenvalues 
exhibit  the  same  near-field  and  far-field  distribution 
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behavior  as  free-space  spatial  eigenvalues.  Spatial  eigen¬ 
value  distributions  for  free-space  and  the  turbulent  atmos¬ 
phere  are  presented  in  the  next  section. 

Eigenvalue  Distributions 

Free-Space  -  The  first  free-sapce,  spatial  eigenvalue 
distribution  presented  is  that  for  a  circularly-apertured , 
coherent  object  field  and  a  circularly-apertured  image  field. 
Circularly-apertured  means  that  the  object  and  image  fields 
are  spatially  limited  to  circular  regions.  Specifically, 
this  means  that  the  object  field  aperture  R.  of  Figure  2  and 
the  image  field  aperture  R 2  of  Figure  2  are  circular  aper¬ 
tures  of  diameters  d1  and  d2  respectively. 

The  second  free-space,  spatial  eigenvalue  distribution 
presented  should  be  that  for  a  circularly-apertured,  incoher¬ 
ent  object  field  and  a  circularly-apertured  image  field.  Un¬ 
fortunately,  this  case  is  not  only  not  found  in  the  literature 
but  it  is  also  difficult  to  derive.  Therefore,  the  second 
free-space  spatial  eigenvalue  distribution  presented  is  that 
for  a si  one-dimensional,  incoherent  object  field  and  a  one¬ 
dimensional  image  field.  One-dimensional  means  that  aper¬ 
tures  R1  and  R2  of  Figure  2  are  slit  apertures  of  widths  d^^ 
and  d2  respectively. 

The  final  free-space  spatial  eigenvalue  distribution 
presented  is  that  for  a  one-dimensional,  coherent  object  field 
and  a  one-dimensional  image  field.  This  type  of  coherent 
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object  field  is  considered  in  addition  to  the  circularly- 
apertured,  coherent  object  field  because  its  spatial  eigen¬ 
value  distribution  is  used  in  the  last  two  sections  of  this 


chapter.  In  all  the  distributions  presented,  the  spatial 
eigenvalues  are  ordered  such  that 

For  the  case  of  a  circularly-apertured,  coherent  object 
field,  Shapiro  has  shown  (Ref  13:2615)  that  the  near-field, 
free-space,  spatial  eigenvalues  exhibit  the  following  rectan¬ 
gular  behavior: 


K-l 


i  ^  D 
i  >  D 


fo 


fo 


i  ^  1  ,  i  £  integers 

(28) 


where 


Dfo  =  (TTd^/AXzY 


(29) 


The  near- field  qualification  of  Eq  (28)  means  that 
>>  1.  For  a  near-field,  circularly-apertured,  coherent 
object  field  (of  diameter  d^)  propagated  a  distance  z  through 
free-space,  D-q  is  called  the  number  of  DOF  of  the  resulting 
circularly-apertured  image  field  (of  diameter  d2).  If  the 
object  field  is  in  the  far-field  (Dfo  <<1),  then  there  is 
only  one  spatial  eigenvalue  and  its  value  is  D^.  In  addi¬ 
tion,  the  image  field  of  a  far-field  object  field  has  only 
one  DOF.  This  would  be  the  case  for  a  point  source. 

Equation  (29)  is  derived  by  combining  Eqs  (26)  and  (28), 
as  is  shown  in  Appendix  C.  Alternatively,  Eq  (29)  can  be 
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derived  using  the  Sampling  Theorem.  The  Sampling  Theorem 
says  (Ref  6:93)  that  the  number  of  DOF  of  an  image  field  is 
obtained  by  dividing  the  solid  angle  subtended  by  the  object 
field  aperture  (when  vie;ved  from  the  image  field  aperture 
R?)  by  the  diffraction  limited  field  of  view  of  the  image 
field  aperture  R  .  The  Sampling  Theorem  derivation  of  Eq  (29) 
is  shown  in  Appendix  D. 

The  spatial  eigenvalue  behavior  given  by  Eq  (28)  is 
illustrated  by  the  distribution  of  Figure  3. 


Figure  3  -  Spatial  Eigenvalue  Distribution  for  Near-Field, 
Circularly-Apertured ,  Coherent  Object  Field 
Propagated  through  Free-Space 

For  the  case  of  a  one-dimensional,  incoherent  object 
field,  Bendinelli,  et.  al.$  have  shown  (Ref  1:1500)  that  the 
near-field,  free-space,  spatial  eigenvalues  exhibit  the 
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following  triangular  behavior: 


X 


<i/Dfo'H 


where 


i  £  Df  '  i  ^  1  ,  i €  integers 


1  >  Dfo' 


(30) 


Dfo'  *  dlVXz 


(31) 


The  near-field  qualification  of  Eq  (30)  means  that 
Dfo'  »  1.  For  a  near-field,  one-dimensional,  incoherent 
object  field  (of  width  d1)  propagated  a  distance  z  through 
free-space,  D~  '  is  called  the  number  of  DOF  of  the  resulting 
one-dimensional  image  field  (of  width  d2).  If  the  object 
field  is  in  the  far  field  (Dfo'  <<  1),  then  there  is  only 
one  spatial  eigenvalue  and  one  DOF  of  the  image  field.  How¬ 
ever,  Bendinelli,  et.  al.,  have  not  shown  that  the  value  of 
the  one  spatial  eigenvalue  is  D^0'» 

The  spatial  eigenvalue  behavior  given  by  Eq  (30)  is 
illustrated  by  the  distribution  of  Figure  4. 

For  the  case  of  a  one-dimensional,  coherent  object  field, 
the  near-field,  free-space,  spatial  eigenvalues  exhibit  the 
following  rectangular  behavior  (Ref  18:801): 

fl  ;  i£Dfo'/2  i^l  ,  i  €  integers 


(  1^0  ;  i  >  Dfo ' / 2  (32) 
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where 


Dfo'/2  =  d1d2/2Xz 


(33) 


Figure  4  -  Spatial  Eigenvalue  Distribution  for  Near-Field, 
One-Dimensional,  Incoherent  Object  Field 
Propagated  through  Free-Space 

The  near-field  qualification  of  Eq  (32)  means  that 
Dfo'/2  >>  1.  For  a  near-field,  one-dimensional,  coherent 
object  field  (of  width  d1)  propagated  a  distance  z  through 
free-space,  D^o'/2  is  called  the  number  of  DOF  of  the  result¬ 
ing  one-dimensional  image  field  (of  width  .  If  the  object 
field  is  in  the  far-field  (Dfo'/2  «  1),  then  there  is  only- 
one  spatial  eigenvalue  and  one  DOF  of  the  image  field. 

The  spatial  eigenvalue  behavior  given  by  Eq  (32)  is  il¬ 
lustrated  by  the  distribution  of  Figure  5. 
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Figure  5  -  Spatial  Eigenvalue  Distribution  for  Near-Field, 
One-Dimensional,  Coherent  Object  Field 
Propagated  through  Free-Space 

Turbulent  Atmosphere  -  A  turbulence  spatial  eigenvalue 
distribution  is  presented  for  a  circularly-apertured,  coher¬ 
ent  object  field  and  a  circularly-apertured  image  field. 
Distributions  are  not  presented  for  a  circularly-apertured 
incoherent  object  field  nor  for  a  one-dimensional  incoherent 
object  field,  as  neither  case  is  found  in  the  literature  nor 
can  they  be  easily  derived. 

For  the  case  of  a  circularly-apertured,  coherent  object 
field,  Shapiro  has  shovm  (Ref  13:2616)  that  the  near-field 
(D^q  »  1),  turbulence  spatial  eigenvalues  are  random  varia¬ 
bles  that  exhibit  with  high  probability  the  following  rectan¬ 
gular  behavior: 


For  a  near-field,  circularly-apertured ,  coherent  object 
field  (of  diameter  d^)  propagated  a  distance  z  through  the 
turbulent  atmosphere,  is  called  the  number  of  DOF  of  the 
resulting  circularly-apertured  image  field  (of  diameter  d^) . 
In  Eq  (35),  X(?2,r^)  represents  the  random  log-amplitude 
fluctuations  of  the  turbulent  atmosphere.  Since  X(?2,r1) 
is  a  random  variable,  both  Df  and  the  number  of  DOF  are  ran¬ 
dom  variables.  If  the  object  field  is  in  the  far-field 
(Dfo  «  1),  then  there  is  only  one  spatial  eigenvalue  and  its 
value  is  D^,.  Equation  (35)  is  derived  by  combining  Eqs  (5C) 
and  (32),  as  is  shown  in  Appendix  E.  When  *<?.?!)  =  0  , 

as  in  free-space,  Eq  (35)  reduces  to  Eq  (29).  The  spatial 
eigenvalue  behavior  given  by  Eq  (34)  is  the  same  as  the  be¬ 
havior  illustrated  in  Figure  3,  but  with  D^q  replaced  by  D^.. 

Energy  conservation  requires  that  the  mean  of  X(r2,r1) 
equal  the  negative  of  the  variance  of  X(r2,ri)  (Ref  4:1374). 
In  Appendix  F  it  is  shown  that  this  requirement  means  that 
■<D^>  =  D^o  .  Since  the  number  of  DOF  of  an  image  field 
is  a  measure  of  the  resolution  capability  of  a  spatial  domain 
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says  that  on 


imaging  system,  the  statement  < 
the  average  the  imaging  system  has  the  same  resolution  capa¬ 
bility  for  the  turbulent  atmosphere  as  it  does  for  free-space. 


Additive  Background  Noise 

This  section  analyzes  the  effect  additive  background 
noise  has  on  the  useable  DOF  of  an  image  field  when  the  ob¬ 
ject  field  is  in  the  near-field,  is  one-dimensional  and  is 
propagated  through  free-space.  In  the  analysis,  apertures 
and  R2  of  Figure  2  are  slit  apertures  of  widths  d^  and 
d2  respectively.  Both  coherent  and  incoherent  object  fields 
are  considered.  The  analysis  is  not  extended,  however,  to 
circularly-apertured  object  fields  nor  to  object  fields  prop¬ 
agated  through  the  turbulent  atmosphere. 

The  analysis  follows  the  example  of  Bendinelli,  et.  al., 
(Ref  1:1499),  and  begins  by  calculating  the  average  IMSE  of 
the  image  field  generated  when  background  noise  is  present 
in  the  image  field  plane  P2  of  Figure  2.  The  background  noise 
is  zero-mean,  spatially  white  noise.  The  average  IMSE  of  the 
image  field  is 


<  IMSE(N)  >  =  /<  [e o(?2)  "  Eo’(?2)] 


>dr„ 


where 


EQ(r2)  = 


(36) 


(21) 
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and  where 


N 

V<V  '  I(bl  +  ni>  <37> 

i=l 

EQ(r2)  is  the  image. field  generated  in  the  absence  of 
background  noise  while  Eo'(r2)  is  an  estimate  of  the  image 
field  generated  when  background  noise  is  present.  The  number 
N  in  Eq  (27)  is  chosen  such  that  the  average  IIISE  of  the  im¬ 
age  field  is  minimized.  This  is  the  reason  N  appears  as  an 
argument  in  Eq  (36).  The  n^  in  Eq  (37)  are  the  random  noise 
coefficients  of  the  zero-mean,  spatially  v/hite  noise.  The 
average  of  Eq  (36)  is  with  respect  to  these  random  noise 
coefficients.  In  Appendix  G  it  is  shown  that  Eq  (36)  is 
equivalent  to 

©o  N 

<  IIISE  (N)  >  =  ^Aihl2  -  J/Xjaj2  -Vn2)  (38) 

i=l  i=l 

? 

where  the  a i  are  object  field  weighting  coefficients  and 
is  the  variance  of  the  background  noise. 

The  next  step  in  the  analysis  is  to  minimize  the  average 
IIISE  given  by  Eq  (38).  In  general,  the  object  field  weight¬ 
ing  coefficients  are  unordered.  Therefore,  in  general  all 
that  can  be  said  is  that  the  second  summation  of  Eq  (38)  must 
be  positive  for  the  average  IMSE  to  be  minimized.  However, 
if  the  object  field  weighting  coefficients  are  ordered  such 
that  ai  ai+i  »  then  the  average  IIISE  is  minimized  by 
choosing  an  N  that  maximizes  the  second  summation  of  Eq  (38). 
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The  second  summation  cannot  be  maximized  by  simply  letting  II 
go  to  infinity  because  of  the  ordering  of  the  spatial  eigen¬ 
values  and  the  ordering  of  the  object  field 

weighting  coefficients  (a^  ^  ai+i^  •  As  N  becomes  too 
large,  | ai |  becomes  less  than  and  the  second  summa¬ 

tion  begins  to  decrease  rather  than  increase.  Therefore, 
the  N  that  maximizes  the  second  summation  of  Eq  (38)  is  the 
N  that  causes  the  term  in  the  summation  to  equal  zero  or  it 
is  the  largest  N  that  insures  the  term  remains  a  positive 
quantity.  The  latter  qualification  is  necessary  because 
the  spatial  eigenvalues  and  object  field  weighting  coeffi¬ 
cients  assume  discrete  values  that  will  not  always  permit  an 
II  to  be  chosen  such  that  the  term  in  the  second  summation  of 
Eq  (38)  is  exactly  zero.  Therefore,  the  N  that  minimizes  the 
average  IMSE  given  by  Eq  (38)  is  the  largest  N  that  satisfies 
the  equation 


Note  that  Eq  (39)  is  equivalent  to  Jbi-nj|2^V^  where 
the  term  |b^_Nj2  is  the  power  of  Nth  spatial  mode  of  the 
image  field  and  is  the  power  (variance)  of  the  additive 
background  noise.  Since  the  N  that  satisfies  Eq  (39)  tells 
how  many  image  field  spatial  modes  have  powers  greater  than 
the  noise  power  (i.e.  these  modes  are  discernible  from  the 
noise),  N  is  the  number  of  useable  DOF  of  the  image  field. 
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Therefore,  Eq  (39)  relates  the  useable  DOF  of  an  image  field 
in  additive  background  noise  to  the  spatial  eigenvalues  of 
the  propagation  channel.  Eq  (39)  is  valid  for  both  coherent 
and  incoherent  object  fields. 

For  a  near-field,  one-dimensional,  coherent  object  field 
propagated  through  free-space,  the  spatial  eigenvalue  be¬ 
havior  is  given  by  Eq  (32)  and  is  illustrated  in  Figure  5. 

The  number  of  DOF  of  the  image  field  when  there  is  no  additive 
background  noise  is  D^q,/2,  which  is  defined  by  Eq  (33).  Be¬ 
cause  of  the  rectangular  behavior  of  the  spatial  eigenvalues, 

D„  1 / 2  is  also  the  N  that  satisfies  Eq  (39).  Therefore,  addi- 
f  o 

tive  background  noise  has  no  effect  on  the  number  of  DOF  of 
the  image  field  when  the  object  field  is  coherent. 

For  a  near-field,  one-dimensional,  incoherent  object  field 
propagated  through  free-space,  the  spatial  eigenvalue  behavior 
is  given  by  Eq  (30)  and  is  illustrated  in  Figure  4.  The  num¬ 
ber  of  DOF  of  the  image  field  when  there  is  no  additive  back¬ 
ground  noise  is  D^0' ,  which  is  defined  by  Eq  (31).  Because 
of  the  triangular  behavior  of  the  spatial  eigenvalues,  Dfo' 
is  not  the  N  that  satisfies  Eq  (39).  The  N  that  satisfies 
Eq  (39)  for  this  case  is 

« ■  Dto'(1  -  ^2/K=n|2>  >  1  (40) 

Equation  (40)  is  obtained  by  letting  i  ■  N  in  Eq  (30) 
and  substituting  the  result  into  Eq  (39).  Since  is 
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the  power  of  the  Nth  spatial  mode  of  the  object  field  and 
2 

Is  the  power  (variance)  of  the  additive  background  noise, 
2  /  2 

the  term  Vrj  /|ai_N|  in  Eq  is  3X1  inverse  power  signal  - 

to  -  noise  ratio  (SNR).  When  there  is  no  additive  background 
noise  (v^  =  0)  ,  the  inverse  SNR  in  Eq  (40)  is  zero  and 

N  =  as  expected.  However,  as  the  additive  background 

noise  increases,  the  number  of  useable  DOF  of  the  image  field 
decreases  until  N  =  0  for  ^  =  ai_N|  .  Therefore, 

additive  background  noise  decreases  the  resolution  capability 
of  an  imaging  system  when  the  object  field  to  be  imaged  is 
incoherent. 

Relationship  Between  Spatial  Domain  and  Frequency  Domain 

In  this  section  it  is  shown  that  when  certain  conditions 
are  satisfied,  it  is  asymptotically  true  that  the  spatial 
eigenvalues  of  a  spatial  domain  imaging  system  are  obtained 
by  sampling  the  power  density  spectrum  (frequency  domain)  of 
the  impulse  response  of  the  imaging  system.  There  are  four 
conditions  that  must  be  satisfied  for  the  above  relationship 
between  the  two  domains  to  be  true.  First,  the  spatial  imag¬ 
ing  kernel  of  the  imaging  system  must  be  space- invariant  so 
that  it  can  be  Fourier  transformed  to  obtain  the  power  density 
spectrum  of  the  system's  impulse  response.  Second,  the  ob¬ 
ject  field  propagated  through  the  system  must  have  a  large 
spatial  extent.  Third,  the  maximum  phase  change  across  the 
object  field  must  be  less  than  some  small  radian  measure  such 
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as  TT/8.  Finally,  the  object  field  must  be  in  the  near-field. 
Near-field  means  the  number  of  DOF  of  the  image  field  asso¬ 
ciated  with  the  object  field  must  be  much  greater  than  one. 

The  relationship  between  the  spatial  domain  and  the 
frequency  domain  is  shown  below  for  a  near-field,  one¬ 
dimensional  object  field  propagated  through  free-space.  Both 
coherent  and  .incoherent  object  fields  are  considered.  For  a 
one-dimensional  object  field,  object  field  aperture  R^  of 
Figure  2  is  a  slit  aperture  of  width  d1 .  For  the  associated 
one-dimensional  image  field,  image  field  aperture  R 2  of 
Figure  2  is  slit  aperture  of  width  d^.  The  propagation  chan¬ 
nel  of  Figure  2  is  free-space. 

For  a  one-dimensional,  coherent  object  field  it  is  shown 
in  Appendix  H  that  the  space-invariant,  coherent  imaging 
kernel  is  of  the  form 

K(  | ?!“?! '  |  )  =  sin(TTd2|r1-r1'| /Xz)/7f  |?1-?1'| 

;  ?1,€R1  (41) 

The  coherent  power  density  spectrum  S  (|?|)  of  the  im- 
pulse  response  associated  with  the  coherent  imaging  kernel 
is  the  Fourier  transform  of  Eq  (41): 


Sc(|?|)  =  rect( |?| A  z/d2) 


(42) 


The  power  density  spectrum  given  by  Eq  (42)  is  illustrated 
in  Figure  6. 


Figure  6  —  Coherent  Power  Density  Spectrum 


The  spatial  eigenvalues  associated  with  the  coherent  ob¬ 
ject  field  are  obtained  by  sampling  Sc(|?|).  If  temporal 
eigenvalues  were  wanted,  the  sampling  frequency  would  be  the 
reciprocal  of  the  time  period  over  which  the  object  field  was 
observed  (Ref  19:206).  Since  spatial  eigenvalues  are  actually 
wanted,  the  sampling  frequency  is  the  reciprocal  of  the  width 
of  the  one-dimensional  object  field.  Since  the  width  of  the 
object  field  is  the  spatial  sampling  frequency  is 
|?x|  =  l/dx  .  Therefore,  for  large  dx  it  is  asymptotically 
true  that  the  spatial  eigenvalues  associated  with  a  near-field, 
one-dimensional,  coherent  object  field  are  given  by 
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K  *  sc(il?ll>  "  so(1/dl> 


i  ^  1,  i€  integers 


(43) 

Applying  Eq  (43)  yields  the  spatial  eigenvalue  distribu¬ 
tion  illustrated  in  Figure  7. 


Figure  7  -  Spatial  Eigenvalue  Distribution  from 

Sampled  Coherent  Power  Density  Spectrum 


The  spatial  eigenvalue  behavior  illustrated  in  Figure  7 
is  in  agreement  with  the  behavior  given  by  Eqs  (32)  and  (33). 

As  the  width  of  the  one-dimensional  object  field  increases, 
the  spatial  sampling  frequency  decreases  and  the  discrete 
distribution  of  Figure  7  approaches  a  continuous  distribution. 

The  above  procedure  is  now  repeated  for  a  one-dimensional, 
incoherent  object  field.  For  this  type  of  field,  Bendinelli, 
et.  al.,  have  shown  (Ref  1:1500)  that  the  space-invariant, 
incoherent  imaging  kernel  is  of  the  form 


Ktjr^r^  )  =  sin2(7Td2  r1~?1 '  | /Xz)/(  7fd2  r1-?1 '  /Xz)2 


;  r1'€R1  (44) 

The  incoherent  power  density  spectrum  Sj(|?|)  of  the 
impulse  response  associated  with  the  incoherent  imaging 
kernel  is  the  Fourier  transform  of  Eq  (44): 

Sjdfl)  =  A-(|?|Xz/d2)  (45) 

The  power  density  spectrum  given  by  Eq  (45)  is  illustrated 
in  Figure  8. 


Figure  8  -  Incoherent  Power  Density  Spectrum 

Again,  the  spatial  sampling  frequency  is  |  J  =  1/dj^ 
Therefore,  for  large  it  is  asymptotically  true  that  the 
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spatial  eigenvalues  associated  with  a  near-field,  one¬ 
dimensional,  incoherent  object  field  are  given  by 


=  S.^ ( i  | T1 1  )  =  Sjti/d^  i  1  ,  i  €  integers  (46) 

Applying  Eq  (46)  yields  the  spatial  eigenvalue  distribU' 
tion  illustrated  in  Figure  9. 


Figure  9  -  Spatial  Eigenvalue  Distribution  from  Sampled 
Incoherent  Power  Density  Spectrum 

The  spatial  eigenvalue  behavior  illustrated  in  Figure  9 
is  in  agreement  with  the  behavior  given  by  Eqs  (30)  and  (31). 
As  the  width  of  the  one-dimensional  object  field  increases, 
the  spatial  sampling  frequency  decreases  and  the  discrete 
distribution  of  Figure  9  approaches  a  continuous  distribu¬ 
tion. 
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The  next  chapter  describes  the  operation  of  two  adaptive 
imaging  receivers  and  compares  the  resolution  capabilities  of 
the  receivers  in  the  spatial  domain  and  the  frequency  domain. 
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IV  Adaptive  Imaging  Receivers 

This  chapter  has  five  purposes.  The  first  purpose  is  to 
describe  the  operation  of  two  ideal,  adaptive  imaging  receiv¬ 
ers.  The  two  receivers  are  the  Channel-Matched  Filter  (CMF) 
Receiver  and  the  Multiplicative-Phase  (MP)  Receiver.  The 
second  purpose  is  to  derive  a  turbulent  coherent  transfer 
function  (CTF)  for  each  receiver.  The  third  purpose  is  to 
derive  a  minimum  average  integrated  mean  square  error  (IMSE) 
expression  for  the  image  field  generated  by  each  receiver. 

The  error  is  the  difference  between  the  image  field  generated 
for  free-space  propagation  versus  the  image  field  generated 
for  propagation  through  the  turbulent  atmosphere.  The  fourth 
purpose  is  to  use  the  derived  turbulent  CTFs  and  minimum  aver¬ 
age  IMSE  expressions  to  compare  the  frequency  domain  and  spa¬ 
tial  domain  resolution  capabilities  of  the  CMF  and  MP  Receiv¬ 
ers.  The  final  purpose  of  this  chapter  is  to  show  the  condi¬ 
tion  that  must  be  satisfied  for  the  results  of  this  chapter 
to  be  applicable  to  object  fields  composed  of  many  isoplanatic 
patches  (extended  object  fields). 

The  first  section  of  this  chapter  presents  the  imaging 
system  geometry  and  field  definitions  used  throughout  the 
chapter.  The  second  section  describes  the  operation  of  the 
CMF  and  MP  Receivers.  The  third  section  derives  a  turbulent 
CTF  and  minimum  average  IMSE  expression  for  the  CMF  Receiver 
while  the  fourth  section  does  the  same  for  the  MP  Receiver. 

The  fifth  section  compares  the  two  receivers  and  the  final 
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section  is  a  discussion  of  extended  object  field  imaging. 

General 

Figure  10  shows  an  adaptive  imaging  system  composed  of 
an  object  field  plane  P1,  receiving  plane  P2  and  an  image 
field  plane  P^.  P^  has  a  circular  aperture  R^  of  diameter 
dl '  P2  ^as  a  c*rcular  aperture  of  diameter  d2 ,  and  P^  has 
a  circular  aperture  R^  of  diameter  d^  (same  diameter  as  R1 ) . 
The  region  between  and  P2  is  called  the  propagation  chan¬ 
nel.  The  region  between  P2  and  P^  is  called  the  adaptive 
imaging  receiver.  The  receiving  plane  P2  is  therefore  both 
the  input  to  the  receiver  and  the  output  of  the  propagation 
channel.  For  mathematical  convenience  the  length  of  the 
adaptive  imaging  receiver  is  the  same  length  as  the  propaga¬ 
tion  channel.  In  a  real  imaging  system  lenses  are  used  to 
considerably  shorten  the  length  of  the  adaptive  imaging  re¬ 
ceiver. 

The  general  impulse  response  of  the  propagation  channel 
of  Figure  10  is  h21(r2,r1).  When  the  propagation  channel  is 
free-space,  the  general  impulse  response  is  replaced  by  the 
free-space  impulse  response  hps^2’^l^'  The  free-sPace  im~ 
pulse  response  is  defined  as 

hpgtrg,?^  =  exp[jk(z  +  |  ?2-r1 1  2/2z)J  /  j  X  z  (47) 
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Equation  (47)  is  based  on  the  paraxial  and  Fresnel  ap¬ 
proximations  (Ref  7:58).  The  paraxial  approximation  is  valid 
since  the  distance  z  between  and  P 2  of  Figure  10  is  much 
greater  than  the  radii  of  R^  and  R2>  The  Fresnel  approxima¬ 
tion  is  valid  since  the  distance  z  between  P^  and  P2  satisfies 
the  following  condition: 


where  (x1,y1)  is  a  point  in  R1  and  (x2,y2)  is  a  point  in  R2 . 


Figure  10  -  Geometry  for  Adaptive  Imaging 


When  the  propagation  channel  of  Figure  10  is  the  turbu- 
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lent  atmosphere,  the  general  impulse  response  is  replaced  by 
the  turbulent  impulse  response  hTB(r2,?1),  The  turbulent 
impulse  response  is  defined  as 


hTB^2' 


x  exp 


[ 


X{r2,r1) 


(49) 


where  X(  ?2,r1)  and  0(r2,r1)  are  turbulence-induced  log- 
amplitude  and  phase  fluctuations  respectively.  They  are 
Gaussian  random  variables  (Ref  17:209)  with  statistics  as 
given  in  Appendix  I. 

EI(r^)  is  defined  as  an  unknown  but  nonrandom  object 
field  in  R1  of  Figure  10.  EQ(r2)  is  deEined  as  the  output 
field  of  the  propagation  channel  of  Figure  10.  By  the 
Huygens-Fresnel  Principle  (Ref  2:Chap  6),  Eo(r2)  is 

Eo(^2)  =  -^EI(?l)h21^2,?l^d?l  (23) 

R1 


EN^2^ 
field  in  R_ 


EN(r2)  is 


is  defined  as  a  zero-mean,  spatially  white  noise 
of  Figure  10.  The  autocorrelation  function  of 


<  en(?2)EN#(?2'  J  >  =  No  ^?2"?2‘  > 


(50) 


Since  EN(r2)  represents  additive  background  noise,  the 
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field  E(r2)  at  the  input  to  the  adaptive  imaging  receiver  is 

e(?2)  =  EQ(r2)  +  EN(r2)  (51) 

The  field  E(r2)  is  vised  in  the  next  section  to  describe 
the  operation  of  the  CMF  and  MP  Receivers. 

Operation  of  Channel-Matched  Filter  and  Multiplicative-Phase 
Receivers 

The  CMF  and  MP  Receivers  are  ideal,  adaptive  imaging 
receivers.  They  are  called  ideal  because  the  general  impulse 
response  h21(r2,?1)  of  a  propagation  channel  is  known  by  the 
receivers.  A  real  receiver  does  not  have  complete  knowledge 
of  this  impulse  response.  They  are  called  adaptive  because 
they  are  able  to  vary  their  response  to  a  propagated  object 
field  such  that  they  partially  compensate  for  the  effects  of 
the  propagation  channel.  The  effects  of  interest  in  this 
thesis  are  turbulence-induced  log-amplitude  and  phase  fluctu¬ 
ations. 

The  CMF  Receiver  operation  is  simply  the  propagation  of 
the  field  E(?2)  from  P2  to  P3  of  Figure  10  using  a  receiver 
impulse  response  that  is  the  conjugate  of  the  impulse  response 
of  the  propagation  channel.  The  resulting  image  field 
ECMF^3^  generated  by  the  CMF  Receiver  is  (Ref  14:2611) 

ECMF^3^  =  J"e(F 2)h2i  (r2,?3)dr2  (52) 
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The  MP  Receiver  operation  is  the  multiplication  of  the 
field  E(?2)  by  an  optimum  phase  and  then  the  propagation  of 
the  field  from  P 2  to  P^  of  Figure  10  using  the  free-space 
impulse  response.  An  optimum  multiplicative  phase  is  a  phase 
that  minimizes  the  average  IMSE  of  the  image  field  generated 

A  _ 

by  the  MP  Receiver.  The  image  field  Ej.jp  (r3)  generated  by  the 
MP  Receiver  is 

EMP(^3)  =  J"E(r2)exp  [j  (j)0(r2)]  hFSCr2,r3)dr2  (53) 

R2  U  ° 

where  0Q(r2)  is  an  optimum  multiplicative  phase  defined  in 
R2  of  Figure  10. 

The  next  section  uses  Eq  (52)  as  the  basis  for  deriving 
the  turbulent  CTF  of  the  CMF  Receiver  and  the  minimum  average 
IMSE  of  the  image  field  generated  by  the  CMF  Receiver. 

Channel-Matched  Filter  Receiver 

The  turbulent  CTF  of  the  CMF  Receiver  is  derived  first 
followed  by  a  derivation  of  the  minimum  average  IMSE  of  the 
image  field  generated  by  the  CMF  Receiver.  All  averages  are 
with  respect  to  the  random  log-amplitude  and  phase  fluctua¬ 
tions  of  the  turbulent  atmosphere. 

Turbulent  Coherent  Transfer  Function  -  When  the  propaga¬ 
tion  channel  of  Figure  10  is  the  turbulent  atmosphere,  the 

*  —  — 

general  impulse  response  h21  (r2,r3)  of  Eq  (52)  is  replaced  ' 

*  _  _ 

by  the  turbulent  impulse  response  hTB  (r2,r3)  defined  by 
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Eq  (49).  For  propagation  through  the  turbulent  atmosphere, 
Shapiro  has  shown  (Ref  14:2612)  that  Eq  (52)  is  rewritten  as 

exp(jk|?3|2/2Z)<ECMF(?3)>  = 

(d2/2\z)/exp(  jk|?1  |2/2z)EI(r1)exp[-D(  | r1-?3 1  > /2] 
x  Jl(7r|V?3|d2/Xz)/(1/|?i-?3|)d?i  (54) 

where  D ( |  J  )  is  the  spherical  wave  structure  function 

(Ref  4:1374).  This  function  is  defined  as 

5 

D(|?l-F3|)  =  1.09k2CN2z(|?1-?3|)3  (55) 

where  k  =  2  7TV  X  is  the  wave  number,  z  is  the  propagation 

2 

distance  of  the  object  field,  and  CN  is  the  refractive  in¬ 
dex  structure  parameter  of  the  turbulent  atmosphere 
(Ref  9:1526).  Equation  (54)  is  derived  in  Appendix  I. 

Shapiro  has  also  shown  that  the  average  image  field 
<ECriF(r3)>  is  obtained  by  convolving  the  object  field 
Ei(?1)  with  the  following  CMF  Receiver  average  impulse  re¬ 
sponse  <hcriF(r)>: 

<hCMF(r)>  =  d2exp[-D(  |r|  )/2  jj1(7T|r|d2/Xz)/2X  z|r| 

(56) 
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Note  that  Eq  (56)  does  not  include  the  quadratic  phase 
factors  of  Eq  (54).  The  factor  on  the  right  side  of  Eq  (54) 
is  not  included  because  the  propagation  distance  z  satisfies 
the  following  condition: 


z  >>  k 


-  I  2 

r, 


MAX 


/  2 


(57) 


Therefore,  the  quadratic  phase  on  the  right  side  of 
Eq  (54)  is  approximately  unity  over  aperture  of  Figure  10 
(Ref  7:61).  The  quadratic  phase  factor  on  the  left  side  of 
Eq  (54)  is  not  included  because  it  is  an  image  field  phase 
that  is  due  to  the  geometry  of  the  imaging  system  rather  than 
to  the  propagation  channel. 

The  average  turbuelnt  CTF  <HCHF(?)>  of  the  CMF  Receiver 
is  obtained  by  taking  the  Fourier  transform  of  Eq  (56).  The 
average  turbulent  CTF  is 


<  H 


CMF 


oo 

(?)>  =  27rJexp[-D(x7rd1/2VD^)/2]j1(27Tx) 


x  J0(27T|H  x)dx 


(58) 


where  Jf’|  =  |f|2Az/d2  and  DfQ  is  defined  by  Eq  (29). 
Equation  (58)  is  derived  in  Appendix  I. 

Shapiro  has  plotted  <HCHF(?)>  for  propagation  through 
the  worst-case  turbulent  atmosphere  (Ref  14:2612).  The 
worst-case  turbulent  atmosphere  means  that  apertures  R1  and 
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R2  of  Figure  10  are  the  same  size  and  the  variance  of  turbu¬ 
lence  log-amplitude  fluctuation  *  .3  0.5.  Shapiro's  plots 
indicate  that  the  average  turbulent  CTF  contains  spatial  fre¬ 
quencies  higher  than  the  free-space,  diffraction-limited  cut¬ 
off  |f'|  .  This  increase  in  spatial  frequency  content  is 
probably  due  to  the  refractive  effect  that  the  turbulent  at¬ 
mosphere  has  ‘on  the  propagated  object  field.  Refraction 
causes  changes  in  the  angle-of-arrival  of  the  propagated  ob¬ 
ject  field  at  the  CMF  Receiver,  which  results  in  a  spreading 
of  the  generated  image  field.  This  spreading  shows  up  as  an 
increase  in  the  spatial  frequency  content  of  the  average  tur¬ 
bulent  CTF.  Because  of  this  effect,  the  spatial  frequency 
cutoff  of  the  CTF  is  not  a  valid  measure  of  the  resolution 
capability  of  the  CMF  Receiver.  In  the  following  subsection 
it  is  shown  that  the  average  IMSE  of  the  image  field  gen¬ 
erated  by  the  CMF  Receiver  is  a  valid  resolution  measure. 

Integrated  Mean  Square  Error  of  Image  Field  -  The  aver¬ 
age  IMSE  of  the  image  field  E^^Cr^)  generated  by  the  CMF 
Receiver  is: 

< IMSE  >  =  /  <  [E(?3)  -  ECMF(r3)]2>d?3  (59) 


where  E(?3)  is  the  image  field  for  free-space  propagation. 
E(?3)  is  obtained  from  Eqs  (23),  (51)  and  (52)  by  replacing 
the  general  Impulse  responses  with  free-space  impulse  re¬ 
sponses.  When  propagation  is  through  the  turbulent  atmosphere 


and  the  additive  background  noise  is  zero-mean,  spatially 
white  noise,  it  is  shown  in  Appendix  J  that  the  average  IMSE 
is 


<  IMSE 


> - /////  E_(u)E_  ( m 

r-»  t»  n  r»  -*•  J* 


_  _  * 

)hFS(v,u)hFS  (n,m) 


x  hps  (v,r3)hpS(n,r3)<  |l-exp  [X(n,m)j 
x  exp  [X(nfr3)J  exp  [-J0<  n,m)  +  j  0(n,r3)] 

+  exp  [X(v,u)  +  X(v,?3  )]  exp  [j0<v.u>] 
x  exp  [-j  0(v  ,r3)]  +  exp[X(v,u)  +  X(n,m)J 
x  exp  [X(n,r3)  +  X(v,r3)]  exp  [j  0(v,u)] 
x  exp  [-j  0(v,r3)]exp[-j  0(n,in)  +  j0(n,r3)]j: 


x  du  dm  dv  dn  dr„  +  N  D_ 

3  o  fo 


;  u,  m€R1  v,  n€R2  (60) 


where  DfQ  is  defined  by  Eq  (29).  The  integral  term  in  Eq  (60) 
is  zero  when  the  object  field  is  isoplanatic.  The  object 
field  is  isoplanatic  when  its  maximum  spatial  extent  is  less 
than  the  atmospheric  coherence  length  rQ  defined  in  Chapter  II. 
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Isoplanatism  permits  the  turbulence-induced  log-amplitude  and 
phase  fluctuations  of  Eq  (60)  to  be  written  as  a  function  of 
only  aperture  R2  of  Figure  10  (Ref  15:463).  For  isoplanatic 
conditions,  it  is  shown  in  Appendix  J  that  Eq  (60)  reduces  to 

<  IMSE  >  =  NoDfo  (61) 

It  is  further  shown  in  Appendix  J  that  Eq  (61)  is  valid 
for  free-space  propagation  even  though  it  is  derived  from 
Eq  (60)  which  is  based  on  propagation  through  the  turbulent 
atmosphere.  Therefore,  when  the  object  field  is  isoplanatic 
and  the  additive  background  noise  is  zero-mean,  spatially 
white  noise,  the  average  IMSE  of  the  image  field  of  the  CMF 
Receiver  is  the  same  for  the  turbulent  atmosphere  as  it  is 
for  free-space. 

Using  a  different  approach  than  that  used  in  Appendix  J , 
Shapiro  has  shown  (Ref  14:2611)  that  the  IMSE  of  the  image 
field,  conditioned  on  knowledge  of  mode-amplitude  estimator 
statistics,  is 

IMSE  =  NQDf  (62) 

where  Df  is  defined  by  Eq  (35).  Since  < Df >  =  Dfo 
(Ref  13:Appendix) ,  Shapiro's  unconditional  IMSE  expression 
is  the  same  as  Eq  (61). 

Shapiro  has  also  shown  that  Eq  (61)  is  not  only  the 
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minimum  average  IMSE  expression  for  the  CMF  Receiver  but  that 
it  is  also  the  smallest  average  IMSE  achievable  by  any  ideal, 
adaptive  imaging  receiver  operating  in  zero-mean,  spatially 
white  background  noise.  The  minimum  average  IMSE  expression 
of  Eq  (61)  is  therefore  a  standard  for  comparing  the  minimum 
average  IMSE  expressions  of  other  ideal,  adaptive  imaging 
receivers. 

The  next  section  derives  the  minimum  average  IMSE  of  the 
image  field  generated  by  the  MP  Receiver  and  the  turbulent 
CTF  of  the  MP  Receiver. 


Multiplicative-Phase  Receiver 

The  minimum  average  IMSE  of  the  image  field  generated 
by  the  MP  Receiver  is  derived  first  followed  by  a  derivation 
of  the  turbulent  CTF„  All  averages  are  with  respect  to  the 
random  log-amplitude  and  phase  fluctuations  of  the  turbulent 
atmosphere. 


Integrated  Mean  Square  Error  of  the  Image  Field  -  The 


average  IMSE  of  the  image  field  EMp(r3)  generated  by  the  MP 


Receiver  is 


< IMSE  > 


(63) 


where  E(?3)  is  the  image  field  for  free-space  propagation. 
E(?3)  is  obtained  from  Eqs  (23),  (51)  and  (53).  In  Eq  (23) 
the  general  impulse  response  is  replaced  by  the  free-space 
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impulse  response  and  in  Eq  (53)  QQ(r^  is  set  equal  to  zero. 

When  propagation  is  through  the  turbulent  atmosphere 
and  the  additive  background  noise  is  zero-mean,  spatially 
white  noise,  it  is  shown  in  Appendix  K  that  the  average  IMSE 
of  the  image  field  is  minimized  if  the  optimum  multiplicative 
phase  0O(?2)  (53)  is  chosen  such  that  the  following 

sufficient  condition  is  satisfied: 


Hike 


*  _  _  _  _ 

)Et  (rnjh^iv  9u)hv^  (n,m)hFS  (v,r3) 


R3R2R1R1 


I'“,iJI  '  7  FS 


x  hFg(n,r3)<  exp  [  X(n,fn)]  exp  [-j  0(n,m)  -  j  0Q ( n) ] 
x[l-  exp  [X(v,u)j  exp  [j  0(v,u)  +  j  0Q(v)j  j  > 


x  du  dm  dv  dr3  >  =  0  >  u,m€R1  v,  n£R2  (64) 


In  Appendix  L  it  is  shown  that,  when  the  sufficient 
condition  above  is  satisfied,  the  minimum  average  IMSE  of 
the  image  field  generated  by  the  MP  Receiver  is 

<  IMSE  >  *  J I J I J  EI(u)EI*(m)hFS(v,u)hFS*(n,in) 


R3R2R2R1R1 


x  hFS  (v,r3)hFg( 


n,r3)£l  -  2<exp  [X(v,u)j 


x  exp  [j0<  V,u)  +  J  0o<v)]>  +  <exp  [X(v,u)j 


/. 

t 


A 


x  exp  [x(n,m)j  exp[j  0(  v,u)  +  j  0Q(v)  -  j  0(n,m)] 
x  exp  [- j  0Q(n)  j  >  j  du  dm  dv  dn  d?3  +  NQDfQ 

;  u,  m  €  R^  v,  n  €  R2  (65) 

Turbulent  Coherent  Transfer  Function  -  Before  the  turbu¬ 
lent  CTF  of  the  MP  Receiver  can  be  derived,  a  choice  must  be 
made  for0Q(r2)  in  Eq  (53).  The  choice  should  satisfy  the 
sufficient  condition  of  Eq  (64)  and  therefore  minimize  the 
average  IMSE  of  the  image  field  generated  by  the  MP  Receiver. 
The  problem  is  that  there  is  no  0Q(r2)  that  meets  this  re¬ 
quirement  without  restrictions  first  being  placed  on  the 
propagation  conditions  or  type  of  object  field  propagated. 
However,  since  many  real  adaptive  receivers  cancel  turbulence- 
induced  phase  fluctuations,  a  compromise  choice  is 
0O<?2)  =  -0<  r2)  .  For  this  choice  of  0Q(r2)  and  Propaga¬ 

tion  through  the  turbulent  atmosphere,  it  is  shown  in  Appen¬ 
dix  M  that  Eq  (53)  is  rewritten  as 

exp(  jkj  r3  j  2/2z)<EMp(r3)  >  =  ( d2/2,X  z)exp(-^j,2/2 ) 

x  ^exp(  jk|r1|  2/2z)EI(?1)exp  [-D^( 0,?1) /2  ] 

Ri 

x  J1(7T|V?3ld2/Xz)(1/l  -?3 1  )d?i  (66) 
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where  D^(0,r^)  is  the  phase  structure  function  defined  by 
Eq  (181)  and  is  the  variance  of  turbulence-induced- 
log-amplitude  fluctuations. 

Since  the  phase  structure  function  of  Eq  (66)  is  not 
space-invariant,  it  is  not  possible  to  obtain  a  space- 
invariant  impulse  response  and  then  Fourier  transform  the 
response  to  obtain  the  turbulent  CTF  of  the  MP  Receiver. 
However,  if  the  object  field  is  isoplanatic,  then  the  phase 

structure  function  in  Eq  (66)  is  zero.  In  this  case,  the 

a  _ 

average  image  field  <  E.Ip(r3)>  is  obtained  by  convolving 
the  object  field  Ej(r1)  with  the  following  IIP  Receiver  aver¬ 
age  impulse  response  <hIIp(r)>: 

<hIIp(r)>  =  exp(-\j^2/2)d2J1(  7fJr|d2/X  z)/2  Xz  |r  |  (67) 

Note  that  Eq  (67)  does  not  include  the  quadratic  phase 
factors  of  Eq  (66).  They  are  not  included  for  the  reasons 
given  in  the  previous  section.  The  average  turbulent 
CTF  <HMp(?)  >  of  the  MP  Receiver  is  obtained  by  taking  the 
Fourier  transform  of  Eq  (67).  The  average  turbulent  CTF  is 

<HMp(f)>  =  exp(-^?r2/2)27rjJ1(27rx)J0(27r|f,|  x)dx  (68) 

o 

where  |?j  -  IN  2  X  z/dg  . 

The  following  section  uses  the  turbulent  CTFs  and 
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minimum  average  IMSE  expressions  derived  in  this  section  and 
the  previous  section  to  compare  the  resolution  capabilities 
of  the  CMF  and  MP  Receivers. 

Comparison  of  Channel-Matched  Filter  and  Multiplicative-Phase 
Receivers 

The  resolution  capabilities  of  the  CMF  and  MP  Receivers 
are  first  compared  in  the  frequency  domain  using  their  turbu¬ 
lent  CTFs.  This  is  followed  by  a  comparison  in  the  spatial 
domain  using  the  minimum  average  IMSE  expressions. 

Since  the  turbulent  CTF  of  the  MP  Receiver  is  based  on 
the  assumption  that  the  object  field  is  isoplanatic,  this 
same  assumption  must  be  applied  to  the  turbulent  CTF  of  the 
CMF  Receiver  before  a  valid  comparison  of  the  two  receivers 
can  be  made.  If  the  object  field  is  isoplanatic,  then  the 
spherical  wave  structure  function  in  Eq  (58)  is  zero  and  the 
average  turbulent  CTF  of  the  CMF  Receiver  is 

<Hcmf<?)>  =  2  7r/j1(2Tr*)J0(2TT|r'|x)dx  (69) 

o 

where  |?'|  =  j  ?  J 2  A  z/d2 

Comparing  Eqs  (68)  and  (69),  it  is  seen  that  except  for 
a  constant  amplitude  attenuation  the  turbulent  CTF  of  the  MP 
Receiver  is  identical  to  that  of  the  CMF  Receiver.  Since 
the  attenuation  in  Eq  (68)  is  independent  of  spatial  frequen¬ 
cy  the  two  turbulent  CTFs  have  the  same  spatial  frequency 
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cutoff.  Since  this  cutoff  is  not  a  valid  measure  of  the  reso¬ 
lution  capability  of  the  CMF  Receiver  it  is  also  not  a  valid 
resolution  measure  for  the  MP  Receiver.  Except  for  the  con¬ 
stant  amplitude  attenuation  in  Eq  (68),  the  turbulent  CTFs 
are  similar  in  another  way .  They  both  approach 
circ(  |  f  |2  Xz/d2)  as  d2~ *■«>,  where  d2  is  the  diameter  of  the 
input  aperture  of  either  the  CMF  or  the  MP  Receiver.  Recall 
that  circ ( | ? } 2 X z/d2 )  is  the  diffraction-limited  CTF  of  a 
free-space  imaging  system  with  a  circular  aperture  of  diame¬ 
ter  d2.  Therefore,  a  CMF  Receiver  with  a  large  input  aper¬ 
ture  achieves  free-space,  diffraction-limited  resolution  of 
isoplanatic  object  fields  propagated  through  the  turbulent 
atmosphere  as  does  a  large-aperture  MP  Receiver  that  only 
cancels  turbulence-induced  phase  fluctuations. 

For  the  joint  conditions  of  propagation  of  an  isoplana¬ 
tic  object  field  through  the  turbulent  atmosphere  and  zero- 
mean,  spatially  white  background  noise,  the  minimum  average 
IMSE  expression  of  the  CMF  Receiver  is  given  by  Eq  (61). 

For  the  same  conditions,  the  comparable  expression  for  the 
MP  Receiver  is  obtained  from  Eq  (65)  by  letting  the  turbu¬ 
lence-induced  log-amplitude  and  phase  fluctuations  be  a  func¬ 
tion  of  only  aperture  R2  of  Figure  10.  This  modification 
of  Eq  (65)  accounts  for  an  isoplanatic  object  field. 

Comparing  Eqs  (61)  and  (65),  it  is  seen  that  the  minimum 
average  IMSE  of  the  image  field  generated  by  the  MP  Receiver  ' 
differs  from  the  minimum  average  IMSE  of  the  image  field 
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generated  by  the  CMF  Receiver.  The  amount  of  difference  is 

4 

the  value  of  the  integral  term  in  Eq  (65).  Since  the  mini¬ 
mum  average  IMSE  of  the  image  field  generated  by  the  CMF 
Receiver  is  the  smallest  average  IMSE  achievable  by  any  ide¬ 
al,  adaptive  imaging  receiver  operating  in  zero-mean,  spa¬ 
tially  white  background  noise,  the  integral  term  in  Eq  (65) 
must  be  non-negative.  Therefore,  the  minimum  average  IMSE 
of  the  image  field  generated  by  the  MP  Receiver  is  equal  to 
or  greater  than  that  of  the  CMF  Receiver.  The  CMF  Receiver 
can  generate  an  image  field  having  degrees  of  freedom 

without  the  average  IMSE  of  the  image  field  exceeding  the 
minimum  average  IMSE  of  Eq  (61).  The  same  is  true  for  the 
MP  Receiver  only  when  the  integral  term  in  Eq  (65)  is  zero. 

There  are  at  least  two  cases  where  a  value  can  be  chosen 
for  0O(?2)  that  not  only  satisfies  the  sufficient  condition 
given  by  Eq  (64)  but  also  eliminates  the  integral  term  in 
Eq  (65).  In  these  two  cases  the  MP  Receiver  generates  an 
image  field  having  the  same  minimum  average  IMSE  as  an  image 
field  generated  by  the  CMF  Receiver.  The  first  case  is 
free-space  propagation  and  the  proper  choice  for  0Q(r2)  is 
0Q(r2)  =  0  .  The  second  case  is  that  of  an  isoplanatic 

object  field  propagated  through  a  turbulent  atmosphere  that 
only  induces  phase  fluctuations  in  the  field.  In  this  case 
the  proper  choice  for  0o(r2)  is  0O(?2)  =  ~  0(?2)  ;  the 

cancellation  of  turbulence-induced  phase  fluctuations.  This 
choice  is  valid  for  all  classes  of  object  fields  (coherent. 
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incoherent,  point  source,  etc.)  as  long  as  the  field  is  iso- 
planatic . 

There  are  at  least  two  cases  where  the  choice  of 
4>0(?2>  -  -4><  ?2)  satisfies  the  sufficient  condition  given 
by  Eq  (64)  but  does  not  eliminate  the  integral  term  in  Eq  (65). 
In  these  cases  the  average  IMSE  of  the  image  field  generated 
by  the  MP  Receiver  is  a  minimum  for  that  receiver  but  it  is 
greater  than  the  minimum  average  IMSE  of  the  image  field  gen¬ 
erated  by  the  CMF  Receiver.  The  first  case  is  propagation 
of  a  point  source  field  through  the  turbulent  atmosphere  while 
the  second  is  propagation  of  an  isoplanatic,  incoherent  object 
field  with  a  real  and  even  object  intensity  distribution  func¬ 
tion  through  the  turbulent  atmosphere.  In  both  cases  the 
turbulent  atmosphere  induces  log-amplitude  and  phase  fluctua¬ 
tions  in  the  propagated  field.  It  is  shown  in  Appendix  N  that 
the  sufficient  condition  given  by  Eq  (64)  is  satisfied  in 
both  cases  when  $Q(r2)  =  “0^2^  *  is  31:50  shoWTl  ln 

Appendix  N  that  the  minimum  average  IMSE  of  the  image  field 
associated  with  a  point  source  field  with  a  complex  ampli¬ 
tude  E„  is 
o 

<  IMSE  >  =  (d2/2  Xz)2EqE0*  -  2exp(-^/2)  +  exp(-^2)j 

*  / [Jl(TT<i2|?3|/Xz)]2(l/|r3|  >2  dr3  ♦  Vfo 
3  (70) 

2 

where  cr  is  the  variance  of  turbulence-induced  log-amplitude 

x 
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fluctuations.  From  Appendix  N,  the  minimum  average  IMSE  of 
the  image  field  associated  with  an  isoplanatic,  incoherent 
object  field  v/ith  a  real  and  even  object  intensity  distribu¬ 
tion  function  I(r^)  is 

<  IMSE  >  =  (d2/2Xz)2[l  -  2exp(-^2/2)  +  exp(-^2) 

I^i)[ji(7rd2|?3-?1|/Xz)]2 

3  1 

x  (1/|?3-?l|)2  d?l  d?3  +  M0Df0  (71) 

The  next  section  presents  the  condition  that  a  imaging 
system  must  satisfy  for  the  results  of  this  section  to  be 
applicable  to  an  object  field  that  extends  spatially  beyond 
an  isoplanatic  patch. 

Imaging  Extended  Object  Fields 

Much  of  the  preceding  work  is  based  on  the  assumption 
that  the  object  fields  do  not  extend  spatially  beyond  a  sin¬ 
gle  isoplantic  patch.  Since  most  object  fields  of  interest 
are  composed  of  many  isoplanatic  patches,  this  would  seem  to 
restrict  the  usefulness  of  the  work.  This  is  not  a  restric¬ 
tion,  however,  provided  the  propagation  channel  of  the  imag¬ 
ing  system  is  underspread,  Shrpiro  has  shown  (Ref  12:472) 
that  an  underspread  channel  is  one  that  insures  that  the  image 
fields  associated  with  spatially  disjoint  isoplanatic  patches 
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of  an  object  field  are  essentially  disjoint  themselves.  Since 
the  image  fields  are  essentially  disjoint,  imaging  an  extended 
object  field  is  reduced  to  an  imaging  of  each  isoplanatic 
patch  and  a  summation  of  all  the  resulting  image  fields. 

The  next  chapter  summarizes  the  conclusions  reached  in 
this  thesis  and  suggests  areas  requiring  further  study. 
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V  Conclusion 


Conclusions 

Space-invariance  in  an  imaging  system  is  a  prerequisite 
for  deriving  frequency  domain  transfer  functions.  In  con¬ 
trast,  the  spatial  domain  counterpart  of  transfer  functions, 
spatial  eigenvalue  distributions,  are  derived  using  the  Normal 
Mode  Approach  without  requiring  an  imaging  system  to  be  space- 
invariant.  This  is  the  main  advantage  of  a  spatial  domain 
imaging  description. 

If  a  spatial  domain  imaging  system  is  space-invariant, 
then  for  a  spatially-large  object  field  it  is  asymptotically 
true  that  the  system's  spatial  eigenvalues  and  hence  its  spa¬ 
tial  eigenvalue  distribution  are  conveniently  obtained  by 
sampling  the  power  density  spectrum  of  the  system's  impulse 
response.  Since  the  power  density  spectrum  is  a  frequency 
domain  concept,  the  two  domains  are  related  through  the  sam¬ 
pling  process.  Besides  space-invariance  and  a  spatially-large 
object  field,  two  other  conditions  must  be  satisfied  for  the 
sampling  to  yield  a  correct  spatial  eigenvalue  distribution. 
The  object  field  being  propagated  through  the  system  must  be 
in  the  near  field  and  the  maximum  phase  change  across  the 
field  must  be  small. 

The  shape  of  an  imaging  system's  spatial  eigenvalue  dis¬ 
tribution  depends  only  on  whether  an  object  field  propagated 
through  the  system  is  coherent  or  incoherent.  If  the  object 
field  is  coherent  the  distribution  has  a  rectangular  shape. 
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If  the  object  field  is  incoherent  the  distribution  has  a  tri¬ 
angular  shape.  The  highest  nonzero  eigenvalue  order  of  a  spa¬ 
tial  eigenvalue  distribution  is  a  function  of  an  imaging  sys¬ 
tem's  input  and  output  apertures,  the  propagation  medium  of 
the  system  ar.d  the  wavelength  of  an  object  field  propagated 
through  the  system.  In  addition,  if  the  object  field  is  in¬ 
coherent  the  highest  nonzero  eigenvalue  order  is  also  a  func¬ 
tion  of  any  additive  background  noise  in  the  imaging  system. 
The  highest  nonzero  eigenvalue  order  is  not  effected  by  ad¬ 
ditive  background  noise  if  the  object  field  is  coherent. 

From  the  frequency  domain,  it  is  concluded  that  a  CMF 
Receiver  with  a  large  input  aperture  achieves  free-space, 
diffraction-limited  resolution  of  isoplanatic  object  fields 
propagated  through  the  turbulent  atmosphere.  The  same  con¬ 
clusion  is  made  for  a  I-IP  Receiver  operating  so  as  to  cancel 
turbulence-induced  phase  fluctuations. 

From  the  spatial  domain,  it  is  concluded  that  for  imag¬ 
ing  isoplanatic  object  fields  in  zero-mean,  spatially  white 
background  noise,  the  CMF  Receiver  earn,  generate  an  image  field 
having  the  maximum  obtainable  DOF  and  also  the  minimum  obtain¬ 
able  average  IMSE  of  any  ideal  adaptive  imaging  receiver. 

There  are  at  least  two  cases  when  the  preceding  statement  also 
applies  to  the  MP  Receiver.  These  cases  are  an  object  field 
propagated  through  free-space  and  an  isoplanatic  object  field 
propagated  through  a  turbulent  atmosphere  that  only  induces 
phase  fluctuations  in  the  field.  In  the  latter  case  the  MP 


Receiver  generates  an  image  field  having  the  minimum  obtain¬ 
able  average  IMSE  by  cancelling  the  turbulence-induced  phase 
fluctuations  in  the  propagated  object  field.  For  point 
source  fields  and  isoplanatic,  incoherent  object  fields  with 
real  and  even  object  intensity  distribution  functions,  phase 
cancellation  minimizes  the  average  IMSE  of  the  image  field 
generated  by  the  MP  Receiver  but  not  to  the  minimum  obtain¬ 
able  average  value  of  the  CMF  Receiver. 

Suggestions  for  Further  Study 

The  relationship  between  the  spatial  domain  and  frequency 
domain  was  shown  for  a  one-dimensional  object  field  propagated 
through  free-space.  This  relationship  should  be  verified  for 
two-dimensional  object  fields  (square  and  circular  apertures) 
and  for  propagation  through  the  turbulent  atmosnhere. 

The  sufficient  condition  given  by  Eq  (64),  for  minimizing 
the  average  IMSE  of  the  field  generated  by  the  MP  Receiver, 
is  complicated.  If  possible,  an  attempt  should  be  made  to 
express  the  optimum  multiplicative  phase  of  the  sufficient 
condition  in  terms  of  the  object  field. 
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Appendix  A  -  Turbulent  Coherent  transfer  Function 

This  appendix  derives  Eq  (13),  which  is  the  average 
transfer  function  for  the  imaging  system  of  Figure  1  when 
the  propagation  channel  is  the  turbulent  atmosphere  and  the 
object  field  is  coherent.  The  derivation  begins  with  Eq  (12): 

<  >  =  ^'0-C0H^<'  exp  [XXXd^fj) 

+  j  0(Xdi|?|)]  >  (12) 

where  %q  COH^^  ^ree“sPace  CTF  defined  by  Eq  (4),  and 

X(Xd^|f|)  and  0(Xd^|?|)  are  respectively  the  log-amplitude 
and  phase  fluctuations  of  the  turbulent  atmosphere.  Since 
X(Xd^|f|)  and  0(Xd^|fj)  are  independent  random  variables 
(Ref  4:1374),  the  average  on  the  right  side  of  Eq  (12)  is 
written  as  two  averages  (Ref  11:211): 

^  ^  =  ^O-COH^^  exp  [  ^di  I  ^1 )] 

x  <  exp[j  0(XdjL|?[)]  >  (1A) 

Prior  to  evaluating  the  averages  on  the  right  side  of 
Eq  (1A),  the  first  and  second  order  statistics  of  X(Ad±|T|) 
and  0(  Adj|  ^  | )  are  given: 

<X(Ad1|f|)  >  -X  -  -r*  (2A) 
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<0(XdJf|)>  =  0 


( 3A) 


<02(Xd.|?j)>  = 


(4A) 


Equation  (2A)  is  a  consequence  of  energy-conservation 
considerations  (Ref  4:1374).  It  can  be  shown  (Ref  4:1374) 
that  the  first  average  on  the  right  side  of  Eq  (1A)  is  equiva¬ 
lent  to 


<  exp  [XCXdJf  |)]  >  =  exp  [<  X(X  di|f|)  >] 

x  expj<  [X(  Xdjf  |  )  -  <  <Y(Xdi|f|  )>]2>/2| 

(5A) 

Substituting  Eq  (2A)  into  Eq  (5A)  yields 

<  exp  [>f<  Xdi|r| )]  >  =  exp(-^2/2)  ( 6A) 

It  can  be  shown  (Ref  4:1374)  that  the  second  average  on 
the  right  side  of  Eq  (1A)  is  equivalent  to 

<  exp  [j  0(  Xdi|?| )]  >  =  exp  [<  j  0(X  djf  | )  >] 

x  exp  j-  <  [0(Xdi|f|)  -  <0(Xdi|?|)>]  >/  2 

( 7A) 

Substituting  Eqs  (3A)  and  (4A)  into  Eq  (7A)  yields 
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<  exp  [X(AaJT|)]>  =  exp(-^2/2) 


(8A) 


Substituting  Eqs  (6A)  and  (8A)  into  Eq  (1A)  yields 


“  ^0-C0H^exp  [—  ^  +V0^2] 


2  2  2 

Finally,,  letting  y"  =  y^  +  yr  yields  Eq  (13): 

<WJ)>  ‘  VcOH(3')eItp'^2/2)  (13) 
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Appendix  3  -  Fredholm  Squation  for  normal  Mode  Imaging 


This  appendix  derives  Eq  (25)  from  Eqs  (19),  (21)  and 
(23).  Eq  (25)  is  then  used  to  obtain  Eq  (26),  which  is  a 
Fredholm  Equation  for  formal  Mode  Imaging.  Equations  (19), 
(21)  and  (23)  are  repeated  below: 


DO 

(19) 

i=l 

oo 

Eo(?2)  =  ^bi'^i^2) 

i=l 

Eo(?2}  =  ZE  I(?l)h21(?2’?l)d?l 

(21) 

(23) 

1 


Substituting  Eqs  (19)  and  (21)  into  Eq  (23)  yields 

oo  r  oo 

^]bi^i(r2)  =  J  ai  0i(r1)h21(r2,r1)dr1 

i=l  Ri  i=l 


(IB) 


Changing  the  order  of  integration-summation  on  the  right 
side  of  Eq  (IB)  yields 


^bi'Fi(r2)  =  $i^i 

i=l  i^l  Ri 


)h21(?2’?l)d?l 


(2B) 


If  the  ^(r^)  and  ^i ( ?2 )  are  chosen  such  that  Eq  (2B) 
is  satisfied  term  by  term,  then  the  following  is  true: 


>i'P1(r2)  =  0i(r1)h21(r2,r1)dr1 


(3B) 
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The  ratio  of  to  a^  is  defined  as 

vai  -  Vxl 


(4B ) 


Substituting  Eq  (48)  into  Eq  (23)  yields  Eq  (25): 

^±^2)  =  J  Q i(?l)h21(?2’?l)d?l 


(25) 


Since  the  ^^(?2)  in  Eq  (25)  are  a  CON  set  of  spatial 
eigenfunctions,  the  following  is  true: 


/'i'i(?2)'i'J*(?2)dr2  = 


(53) 


Expressing  Eq  (25)  in  terms  of  ^(r^)  and  the  conjugate 
of  Eq  (25)  in  terms  of  (?2),  and  then  substituting  the 

resulting  expressions  into  Eq  (58)  yields 


)h21(?2’?l)d?l 


-  <WMj  >// <?!<?! 

x  /^j‘(?l’)h2i*(?2'?l')d?l  d?2  !  ?l'eRl 

< 

Changing  the  order  of  integration  in  Eq  (6B)  yields 

-  (1/V  Aj)/ //hnt?2>?l)h2i*(?2’?l')d?2 


R1R1R2 


67 


68 


Appendix  C  -  Decrees  of  Freedom  Expression  for  a 


Near-Fidd,  Circularly-Apertured ,  Coherent 
Object  Field  Propagated  through  Free-Space 

This  appendix  derives  Eq  (29),  which  is  the  DOF  expression 
for  a  near-field,  circularly-apertured ,  coherent  object  field 
propagated  through  free-space.  The  derivation  begins  with 
Eq  (25): 


Ai*i(V)  =  I  K(r1',r1)  0i(?1)dr1  ;  '  €  R1 

R. 


(25) 


Multiplying  both  sides  of  Eq  (25)  by  0i*(r1')  and  inte¬ 


grating  over  R1  yields 


R-  R1R1 


x  d?1  dr^  ;  ?1€R1  (1C) 


Using  the  orthogonality  of  the  (?1)  to  simplify  the 
left  side  of  Eq  (1C)  yields 


\  =  ffK(  r1' ,rx)  0±(r1)  01*(r1<  )dr1  d?1  *  ;  €  R 

R1R1 


(2C) 


Summing  both  sides  of  Eq  (2C)  over  all  values  of  i  yields 

=  //*<? X0i(?i>  0i#(?i’)d?i  drx' 


R1P! 


1=1 
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;  rl  €  R1 


(3C) 


Applying  Mercer's  Theorem  to  the  right  side  of  Eq  (3C) 


yields 


IK  =  I J  '  ,r^)  J^r^'-r^dr^  d?x  •  ;  r1 '  € 


i=l  R1R1 


Using  the  sifting  property  of  the  delta  function,  Eq  (4C) 
is  simplified  to 


XjV  =  J  K(r1,r1)dr1 


For  free-space  propagation,  the  imaging  kernel  KCr^,?^) 


in  Eq  (5C)  is  defined  as 


K( ?i , ?! )  = 


Hpg  ( r g  ^ 2_  ^  ^ps  ^  ^ 2  *  ^  1  ^  ^^2 


where 


hpgtrg.i^)  =  exp[jk(z  +  |  Fg-rJ  2/2z)]/ j  Xz 
Substituting  Eq  (47)  into  Eq  (6C)  yields 
K(?1,r1)  =  (l/Xz)2/dr2  =  Tf(d2/2Xz)2 


where  R2  is  a  circular  aperture  of  diameter  dg.  Substituting 
Eq  (7C)  into  Eq  (5C)  yields 
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(8C) 
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Appendix  D  -  Sampling  Theorem  Derivation 


of  a  Degrees  of  Freedom  Expression 

This  appendix  derives  Eq  (29)  using  the  Sampling  Theorem. 
Equation  (29)  is  the  DOF  expression  for  a  near-field,  circu- 
larly-apertured,  coherent  object  field  propagated  through 
free-space.  The  Sampling  Theorem  says  (Ref  6:93)  that  the 
number  of  DOF  of  am  image  field  is  obtained  by  dividing  the 
solid  angle  subtended  by  an  object  field  aperture  R^  (when 
viewed  from  the  image  field  aperture  R2)  by  the  diffraction 
limited  field  of  view  of  the  image  field  R2«  In  this  appen¬ 
dix  apertures  R1  and  R2  are  circular  apertures  of  diameters 
dj^  and  d2  respectively.  The  apertures,  solid  angle  and 
field  of  view  0DL  mentioned  above  are  illustrated  in  Figure  11. 

of  Figure  11  is  an  object  field  plane  while  P2  is  an  image 
field  plane.  The  planes  are  separated  by  free-space. 

The  solid  angle  subtended  by  the  object  field  aper¬ 
ture  R^  of  Figure  11  is 

Qs  =  A^/z2  (ID) 

where 

ARl  =  7Td12/4  (2D) 

The  diffraction  limited  field  of  view  QDL  of  the  image 
field  aperture  Rg  of  Figure  11  is  (Ref  6:93) 
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* 


(3D) 


n 


DL  = 


where 


ar2  -  7Td22/4 


(4D) 


Figure  11  -  Sampling  Theorem  Geometry- 


Dividing  Eq  (ID)  by  Eq  (3D)  yields  Eq  (29): 

Dfo  “  (ARi/z2)/(  X2Mr2)  =  (Ttd^/AXz)2  (29) 

where  D is  the  number  of  DOF  of  the  image  field  associated 
fo 

with  the  object  field  described  at  the  beginning  of  this 
appendix. 
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Appendix  £  -  Degrees  of  Freedom  Expression  for  a 
Near-Field.  Circularly-Apertured,  Coherent  Object 
Field  Propagated  through  The  Turbulent  Atmosphere 


This  appendix  derives  Eq  (35),  which  is  the  DOF  expres¬ 
sion  for  a  near-field,  circularly-apertured,  coherent  object 
field  propagated  through  the  turbulent  atmosphere.  The  deri¬ 
vation  begins  with  Eq  (5C)  from  Appendix  C: 


,r*  )dr1 


( 5C  ) 


For  propagation  through  the  turbulent  atmosphere,  the 
imaging  kernel  K(?1,r1)  in  Eq  (5C)  is  defined  as 


K< 


^1*^1^  ~  j *^l^hTB  ^2'^l^d^2 


where 

hTB^2*^l^  =  jexp[jk(z  +  |  ?2-r  J  2/2z)]  /  j  X  z| 


(IE) 


x  exp[X(i'2,r1)  +  j0(r2,?1)]  (49) 


Substituting  Eq  (49)  into  Eq  (IE)  yields 


K( 


rlf?i)  =  (1/ Xz)2/exp[2  Xd^r^Jdr, 


(2E) 


Substituting  Eq  (2E)  into  Eq  (5C)  yields 
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=  (l/\z)2f  f  exp[2  >^(r2,r1)]  d?2  drx 
i=i  RiR2 


(3E) 


For  a  near-field,  circularly-apertured,  coherent  object 
field  propagated  through  the  turbulent  atmosphere,  the  spa¬ 
tial  eigenvalue  behavior  is  given  by  Eq  (34): 


v  _,1  i  i<.Dt 

''i  “ 

0  ;  i  >  D, 


i  ^  1  ,  i  €  integers 


(34) 


From  Eq  (34),  the  summation  of  Xi  over  all  values  of  i 
yields 


oo 


2X  -  d, 


(4E) 


i=l 


Equating  Eqs  (3E)  and  (4E)  yields  Eq  (35): 
Df  =  (l/Xz)2JP-fexp[2X(r2,r1)]dr2  d?1 


(35) 


R1R2 
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Appendix  F  -  The  Consequence  of  Energy  Conservation  with 


Regard  to  the  Degrees  of  Freedom  of  an  Image  Field 

This  appendix  shows  that  since  the  mean  of  Xlr  ,rx) 
must  equal  the  negative  of  the  variance  of  X(  r2,?1),  the 
average  of  Df  must  equal  D^.  X(i'2,r1)  represents  the  log- 

amplitude  fluctuations  of  the  turbulent  atmosphere.  is 
the  number  of  DOF  of  the  image  field  associated  with  a  near- 
field,  circularly-apertured,  coherent  object  field  propagated 
through  the  turbulent  atmosphere  while  Df  is  the  number  of 
DOF  of  the  image  field  associated  with  the  same  object  field 
propagated  through  free-space.  <  Df  >  and  Dfo  are  defined 
as  follows: 


Df  =  (l/\z)2J  J  <  exp  [2  X(r2,r1)]>dr2  dr1  (IF) 
R1R2 

Dfo  =  (1/A  z)2//  d?2  dr^  (2F) 


R1R2 


It  can  be  shown  (Ref  3:71)  that  the  average  in  Eq  (IF) 
is  equivalent  to 


exp  [2  X(?2,?1)]>  =  exp[2<X(?2,r1)>] 

x  exp|<[2X(?2,r1)  -  2  <  X(r2,?1)  >  ]  2  >/2  \  (3F) 


Energy  conservation  requires  the  following  to  be  true 


(Ref  4:1374): 


<X(?2,?1)>  •]  =  -^r2  (4F) 

where  £  is  the  mean  of  and  is  the  variance  of 

X(r2,?i).  Substituting  Eq  (4F)  into  Eq  (3F)  yields 

<exp[2X(?2,r1)j  >  =  1  (5F) 

Substituting  Eq  (5F)  into  Eq  (IF)  reduces  Eq  (IF)  to 
Eq  (2F)  and  yields 

<Df  >  *  Df0  (6F) 
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;e  Field  in  Additive  Background  Noise 


This  appendix  derives  Eq  (38),  which  is  the  average  IMSE 
of  a  Normal  Mode  image  field  in  additive  background  noise. 

The  derivation  follows  the  example  of  Bendinelli,  et.  al.t 
(Ref  1:1499).  The  average  IMSE  of  the  Normal  Mode  image 
field  is 


< IMSE(N) > 


-  /<[Eo<V  - 


where 


w  - 

i=i 

and  where 


1) 

£(b1»n1)'P1(?2) 


(21) 


(37) 


EQ(r2)  is  the  image  field  generated  in  aperture  R2  of 
Figure  2  in  the  absence  of  background  noise  while  E0'(r2)  is 
an  estimate  of  the  image  field  generated  when  background 
noise  is  present.  The  number  N  in  Eq  (37)  is  chosen  such 
that  the  average  IMSE  of  the  image  field  is  minimized.  This 
is  the  reason  N  appears  as  an  argument  in  Eq  (36).  The  n^ 
in  Eq  (37)  are  the  random  noise  coefficients  of  the  zero-mean, 
spatially  white  noise.  The  average  of  Eq  (36)  is  with  respect 
to  these  random  noise  coefficients.  The  ¥i(r2)  in  Eqs  (21) 


and  (37)  are  a  CON  set  of  spatial  eigenfunctions  defined  in 


image  field  aperture  R2  of  Figure  2.  The  in  Eqs  (21)  and 
(23)  are  image  field  weighting  coefficients  that  are  related 
to  object  field  weighting  coefficients  by 


>i  =  VaT  ai 


(1G) 


Substituting  Eq  (1G)  into  Eqs  (21)  and  (37),  and  then 
substituting  Eqs  (21)  and  (37)  into  Eq  (36)  yields 

< IMSE(N) >  =  /<  F,<V 

R2  L i=N+l 


N  *i2 

i  =  l 


(2G) 


Expanding  Eq  (2G)  and  reversing  the  integration-summation 
order  yields 

<IMSE(N)>  =  £ 

i=N+l  j=N+l  2 

+  ^  ^  <ninJ  >  i(V  'J'/(?2,dr2 

i=l  j=l  2 

-  2Z  ^V\ai<nJ>  /^i<?2,^j*<?2)d?2 


issN+1  j=l 


(3G) 


By  the  orthogonality  of  the  ) 


'ij 
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(4G) 


Substituting  Ea  (4G)  into  Eq  (3G).  and  noting  that  the 
last  term  in  Eq  (3g)  is  zero  because  the  summations  are  over 
disioint  intervals,  yields 


< IMSE(N) > 


oo 


Y.  K 


i=N+l 


N 


2  t 

— » 

"i  +  ] 

_  <ni 

i=l 


(5G) 


Since  the  additive  background  noise  is  zero-mean,  the 
average  on  the  right  side  of  Eq  (5G)  is 


<  |  ni 


2^  _  ^2 


(6G) 


o 

where  rr-  is  the  variance  of  the  noise.  Substituting  Eq  (6G) 
into  Eq  (5G)  yields  (after  a  rearrangement  of  terms)  Eq  (38): 

oo  N 

<  IHSE(N)  >  =  ^<Aih|2-^)  (38) 

i=N+l  i=l 
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One-Dimensional  Coherent  Object  Field 
Propagated  through  Free-Space 


This  appendix  derives  Eq  (41),  which  is  the  space- 
invariant  imaging  kernel  for  a  near-field  (D^q'/2  >>  1), 
one-dimensional,  coherent  object  field  propagated  through 
free-space.  With  respect  to  Figure  2,  the  general  expression 
for  the  imaging  kernel  K(r1*,r1)  is  given  by  Eq  (27): 


K( 


iV  .Tf)  =  / 


*  —  — 


21 ( ^2  *^1 )h2l  ^r2,?l‘^dr2  ’  ^l'€Rl  (27) 


where  is  the  object  field  aperture, R2  is  the  image  field 
aperture  and  h21(?2,r1)  is  the  general  impulse  response  of 
the  propagation  channel.  For  a  one-dimensional  object  field, 
is  a  slit  aperture  of  width  d^.  For  the  associated  one¬ 
dimensional  image  field,  R2  is  a  slit  aperture  of  width  d2« 
For  free-space  propagation,  the  general  impulse  response 
h2i(r2,?i)  is  replaced  by  the  free-space  impulse  response 
hFS^2’^l)*  For  a  sPace_invarian‘t  imaging  kernel,  K(?1',r1) 
is  replaced  by  K(|r1-r1,|  ).  Therefore,  for  a  near-field 
(DfQ,/2  »  1),  one-dimensional,  coherent  object  field  propa¬ 
gated  through  free-space,  the  space -invariant  imaging  kernel 
from  Eq  (27)  is 


where 


hFS(?2,r1)  =  exp[jk(z  +  |  |  2/2z)] /j  X  z  (47) 

Substituting  Eq  (47)  into  Eq  (1H)  yields 

K(  |  '  |  )  =  (1/Xz)2exp  [jk(|?1|2  -  |r1'j2)/2zj 

+d2/2 

x  J  exp[-jkr2.(r1-?1,)/2zjdr2  ; 

“d2^2  ( 2H ) 

Evaluating  the  integral  in  Eq  (2H)  and  using  the  iden¬ 
tity  sin(w)  =  £exp(w)  -  exp(-w)J/2j  yields 

K(|?l-?i»|)  =  (1/Xz)exp[jk(|r1|2  -  |?1,|2)/2z] 
x  sin(7Tcl2|r1-?1,| /Xz)/7f 

;  ^i'eRi  (3H) 

If  the  maximum  phase  change  across  the  object  field  is 
less  than  some  small  radian  measure  such  as  7T/8,  then  the 
quadratic  phase  term  in  Eq  (3H)  may  be  eliminated.  The  condi¬ 
tion  for  elimination  is 

K(l?i|2  -  I?i'|2,Ihax/2z  «  ma  (4H) 
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Since  the  width  of  the  object  field  aperture  is  d^,  the 
quantity  ([r^j2  -  |  *  |  2 )  j  max  in  Eq  is  ecluivalent  to 

^lrl|  “  |rl'|  ^|mAX  =  (di/2^2  (5H) 

Substituting  Eq  (5H)  into  Eq  (4H)  yields 

z  »  2d*/\  (6H) 

When  Eq  (6H)  is  satisfied,  Eq  (3H)  reduces  to  Eq  (41): 
K(  | ?!-?!  •  |  )  =  sin(7Td2|r1-r1*| /Xz)/7f  jr^rj^'l 

;  r1'€R1  (41) 

where  the  (1/Xz)  constant  in  Eq  (3H)  is  dropped  because  it 

is  unimportant  to  the  general  form  of  the  coherent  imaging 
kernel. 

In  addition  to  Eq  (6H),  there  is  another  condition  that 
must  be  satisfied  for  Eq  (41)  to  be  true.  The  additional 
condition  is  that  the  one-dimensional  coherent  object  field 
be  in  the  near-field.  This  means  that 

Dfo'/2  ^  1  (7H) 
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where 


Dfo'/2  =  d1d2/2\z  (33) 

Substituting  Eq  (33)  into  Eq  (7H)  yields 

d2  »  2  Xz/d1  (8H) 

Therefore,  for  Eq  (41)  to  be  true,  the  propagation  dis¬ 
tance  z  must  be  large  enough  to  satisfy  Eq  (6H)  and  the  width 
of  the  image  field  aperture  must  be  large  enough  to  satisfy 
Eq  (8H). 


( 
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Appendix  I  -  Image  Field  and  Average  Coherent  Transfer 


Funct i on  of  the  Channel-Hatched  Filter  Receiver 


This  appendix  derives  Eqs  (54)  and  (58).  Equation  (54) 
is  the  expression  for  the  image  field  generated  by  the  CMF 
Receiver  while  Eq  (58)  is  the  expression  for  the  average  CTF 
of  the  CMF  Receiver.  The  derivation  of  both  equations  begins 
with  Eq  (52): 


ECMF(r3^ 


=  /  E(r 


2)h2i*(?2,r3)d?2 


If  propagation  is  through  the  turbulent  atmosphere  then 
the  general  impulse  response,  h21(r2,r3),  of  Eq  (52)  is  re¬ 
placed  by  the  impulse  response,  hTB(r2,?3).  Equation  (52) 
is  then  rewritten  as 


ECMF^3^  -  j E^2^hTB  ^2,r3^d^2 
R2 


where 


hTB(r2,r3)  =  |exp[jk(z  +  J  ?2~?3|  2/2z)] /J  X 


x  exp  [X(r2,r3)  +  j0(r2,r3)] 


The  field  E(?2)  of  Eq  (52)  is  the  input  to  the  CMF 
Receiver.  By  Eq  (51)  it  is 


< 


(51) 


e(?2)  =  eq(?2)  +  en(?2) 


where  EQ(r2)  'the  outPut  field  of  the  propagation  channel 
and  En(?2)  is  a  zero-mean,  spatially  white  noise  field  repre¬ 
senting  additive  background  noise.  For  propagation  through 
the  turbulent  atmosphere  the  field  EQ(?2)  is  expressed  as 


Eo(?2> 


-I 

R, 


El(?l)hTB(?2,?i)d?i 


(21) 


where  EjCr^)  is  an  unknown  but  nonrandom  object  field  in 
aperture  R1  of  Figure  10.  Substituting  Eqs  (51)  and  (21) 
into  Eq  (II)  yields 


A 

E, 


MF ^ r3 ^  =  ^ -^RI^rl^TB^r2*rl^TB  ^r2,r3^  ^rl  ^r2 
R2R1 

/en<7; 


2)hTB  (?2,r3)dr2 


(31) 


Since  X(r',r)  and  0(r,,r)  are  random  variables,  the 
turbulent  impulse  response  given  by  Eq  (49)  is  a  random  func¬ 
tion.  Therefore,  the  image  field  given  by  Eq  (31)  is  also 
random.  On  the  average  however  the  field  is 


A 

<  E 


CMF^r3^ 


.// 

R2R1 


E j ( r^ ) ^ h ( r2 , r^ ) h (r2,r3)^^ 


x  dr^  d?2 


(41) 
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where  the  second  term  in  Eq  (31)  is  dropped  because  the  noise 
field  is  zero-mean.  Substituting  Eq  (49)  into  Eq  (41)  yields 


exp(  jk|r3|2/2z)<ECMF(r3)>  =  ( 1/X  z) 2J"exp(  jk | r1 1  2/2z) 

R1 

'5^)  J exp  ^-jk(?1-r3)  ,?2/zj<  exp  [*(?.,.?!>] 


X  EI(: 


x  exp  [X(r2,r3)]  ><exp  [j  0  (^2,?^^)  -  j0(?2,r3)]> 


x  dr2  dr1 


(51) 


where  the  average  of  the  right  side  of  Eq  (51)  is  written  as 
two  separate  averages  of  X(r',r)  and  0(r',r)  since  these 
are  independent  random  variables  (Ref  4:1374).  Furthermore, 
X(r',r)  and  0(r',r)  are  Gaussian  random  variables 
(Ref  4:1374)  with  the  following  statistics: 


<X(r',r)>  =  £  =  - 

(61) 

<  0( r*  ,r)  >  =  0 

(71) 

<02(r'  ,r)>  = 

(81) 

It  can  be  shown  (Ref  4:1374)  that  the  first  average  on 
the  right  side  of  Eq  (51)  can  be  expressed  as 


( 


87 


<cxp[X(r2,rl)  +  X(?2.r3)]>  =  exp  [<  X(?2  tr± )  >] 

* 

x  exp  [<  X(?2  f  ?3)>]  exp<  <  [X( ?2  ,r1) 

* 

"  <X(r2,r1)>]2>/2>expj<[X(r2,r3) 

-  <X(?2,?3)>]2>/2| 


Substituting  Eq  (61)  into  Eq  (91)  yields 
<exp  [XXr^tr^)  +  X(r2,r3)]>  = 

exp(-2^)exp|[<X2(r2,r1)>-  2]/2^ 

x  exp|[<X2(?2,?3)>  -  X2]/^ 


(101) 


Since  =  <<X'2(F',F)>  - ,  Eq  (101)  is  simplified 


<  exp  +  X(r2,r3)]>  = 


exp<-[<X2(?2»?l^>  +  <X2(?2,?3)>  -  2i2 


Equation  (111)  is  reduced  further  to 


« 

]/2  i 

(111) 


.?*?«<*> — ; 
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(121) 


:  ( 


<  exp  [Xt^,?^  +  X(?2,r3)]>  = 


>  * 

"  -  <  [X^,?^  -  X(r2,?3)J2>/2  ► 
> 


Now,  the  log-amplitude  structure  function  is  defined  as 


(Ref  4:1374) 


Dx(p',p)  =  <[X(r*+p‘  ,  r+p)  -  X(r',r)J2>  (131) 


Using  Eq  (131)  in  Eq  (121)  where  p'  =0  and 
p  ■  ?1-r3  gives 


<exp  [X(r2,r1)  +  X(?2,r3)]>  =  exp[-Dx(  J  r.j-r3|  ) /2  J 


(141) 


The  second  average  on  the  right  side  of  Eq  (51)  can  be 


expressed  as 


<  exp[j  0(r2,r1)  -  j0(?2,r3)]>  = 

exp  [  J  <0(?2,r1)>  -  j  <  0(r2,r3)>] 
x  expj-  <[0(?2,r1)  -  <  0(r2tr1)  >j2>/s| 
x  expj-  <  [0(?2,r3)  -  <  0(?2,r3)  >j2>/2  ►  (151) 


Substituting  Eq  (71)  into  Eq  (151)  gives 


<exp  [j  -  j$(r2,r3)]>  = 

expj- [<  02(?2,r1)>  +  <02(r2,?3)>]/2 


>  (161) 


Since  0(r',r)  is  zero-mean,  Eq  (161)  is  simplified  to 
<exp[j  (f)(r2,r1)  -  j0(r2,r3)]>  = 

<-<[0(?2,?l)  -  0(?2,?3)]2>/2| 


exp< 


(171) 


Now,  the  phase  structure  function  is  defined  as 
(Ref  4:1374) 


D^P'.P)  =  <[0(r,+pl  ,  r+p)  -  0(r’,r)]2>  (181) 

Using  Eq  (181)  in  Eq  (171)  where  p'  =  0  and 
p  =  ?1~r3  yields 

< exp  [j  0(r2,rx)  -  j  0(?2,r3)]>  »  exp  [-D^(  l^-r^  )/2] 

(191) 

Substituting  Eqs  (141)  and  (191)  into  Eq  (51)  gives 
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t 


i 


exp(  jk|r3J2/2z)<ECHF(r3)>  = 

(1/X  z)2J"exp( jk|r1| 2/2z)Ej( ?1 )exp [-D( J r.j-r3 j )  / 2^ 
R1 

x  J*exp[-j2  7T(r1-r3)  .r2/Xz]dr2  d?1 


(201) 


where  D(|rirr3|)  =  D^l^-r^)  +  D0( | )  is  the 
spherical  wave  structure  function  and  where  the  wave  number 
k  in  the  inner  integral  is  replaced  by  2  7T/X.  Note  that  the 
inner  integral  of  Eq  (201)  is  the  two-dimensional  spatial 
Fourier  transform  of  aperture  R2 ,  where  (r^-?3  )  /X  z  —  f 

is  a  two-dimensional  spatial  frequency  vector.  Since  R2  is 
a  circular  aperture  of  radius  d2/2,  the  inner  integral  is 
actually  a  Fourier-Bessel  transform  of  circ(2 j?2 | /d2)  where 

1  !  l?2|^V2 


circ( 2  r2  /d2)  =  < 


otherwise 


(211) 


From  a  table  of  Fourier-Bessel  transforms,  the  transform 
of  Eq  (211)  is 

Jcirc  ( 2  |  ?2  |  /d2 )  j  =  d2XzJ1(7T|r1-r3|d2/Xz)/2|r1-r3j 

(221) 

Replacing  the  inner  integral  of  Eq  (201)  with  Eq  (221) 
yields  Eq  (54): 


( 
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? 


I 


exp(jk|r3|2/2z)<EcriF(r  )>  = 


(d2/2  Xz)J"exp(  jk  |r1 | 2/2z)EI(r1 )exp[-D(  |  r1-?3|  )/2j 
R1 

x  Ji(  Tf|r1-r3|.d2/X  z)/(l/|  Fj-Fg  J  )d?x  (54) 


If  the  quadratic  phase  factors  in  Eq  (55)  are  ignored, 

A  _ 

then  the  CMF  Receiver  average  image  field,  <  ECMF^r3^ -->*  is 
obtained  by  convolving  the  object  field,  E-j-tr^),  with  the 
following  CMF  Receiver  average  impulse  response  <hCIIp(r)>  : 

<ChcMF(*;)>  =  d2exp[-D(  |r|)/2]j1(TT|r|d2/Xz)/2Xz  r 

(56) 

The  average  CTF,  «C  HCMF(? )  >,  of  the  CMF  Receiver  is 
obtained  by  taking  the  Fourier  transform  of  Eq  (56): 

+  OO 

<hCMf^)>  =  (d2/2\z)  J'exp[-D(  jr|  )/2]jx(  7T|r|d2/X  z) 

—  oo 

x  exp(-j2  Tf|f  |  |r j ) ( 1 / 1 r | )dr  (231) 

Let  x  =  |r|d2/2Xz  and  |  ?'  j  =  |?|2Xz/d2  so  that 
Eq  (231)  becomes 

+oo 

<HCMF^^  =  J*exp[-D(x2Xz/d2)/2  jj1(2  7Tx) 

~0O 

x  exp(-j2  TTj?' J  x) ( l/x)dx  (241) 
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Expressing  the  argument  of  the  spherical  wave  structure 
function  in  terms  of  DfQ,  the  free-space  DOF  for  a  circularly- 
apertured  object,  yields 

+  oo 

<HCMF(?)>  =  (d2/2^z)  Jexv[-D(xTCd1/2'^D^))/2]j1(2TTx) 

—  oo 

x  exp(-j27T|f *|  x)(l/x)dx  (251) 

where 

Dfo  =  (7Td1d2/4Xz)2  (29) 

Using  a  table  of  Fourier-Bessel  transforms  (Ref  11:145), 
Eq  (251)  is  simplified  to  Eq  (58): 

f  OO 

<HCMF(f)>  -  2  7T Jexp[-D(xTrd1/2YD^)/2]j1(27Tx) 

o 

x  Jq(2  7T|?' j x)dx  (58) 
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Appendix  J  -  Average  IMSE  of  the  Image  Field  of  the 
Channel-Matched  Filter  Receiver 

This  appendix  derives  the  expression  for  the  average 
IMSE  of  the  image  field  generated  by  the  CMF  Receiver,  when 
propagation  is  through  the  turbulent  atmosphere  and  the  addi¬ 
tive  background  noise  is  zero-mean,  spatially  white  noise. 
Another  average  IMSE  expression  is  derived  for  the  same  con¬ 
ditions  in  conjunction  with  isoplanatism.  For  the  background 
noise  described  above,  it  is  shown  that  the  average  IMSE  for 
the  joint  conditions  of  isoplanatism  and  propagation  through 
the  turbulent  atmosphere  is  the  same  as  the  average  IMSE  for 
free-space  propagation.  The  average  IMSE  of  the  image  field 
generated  by  the  CMF  Receiver  is 

< IMSE  >  =  /<  [e(?3)  -  W?3)]  2>d?3  (59) 

R3 


where  E(?3)  is  the  field  for  free-space  propagation  and  is 
expressed  as 


*  //e I<?1 

Vi 

*  /en(? 


E^r3^  =  J  J  EI^rl^hFS^2,*'l^hFS  ^2’^3Jcrl  dr2 
R2R1 


*  _  _  _ 

2)hFS  (r2,r3)dr2 


(1J) 


Using  Eqs  (23),  (51)  and  (52),  the  image  field  ECMF(r3) 


generated  by  the  CMF  Receiver  is  expressed  as 


ECMF^r3^  =  f  jE I^rl^h21^r2»rl^h21  ^2’^Vdrl  d^2 
R2R1 


+  lE  N<? 

R„ 


NVi2^h21  ^r2’^3^d^2 


(2J) 


In  Eqs  (1J)  and  (2J),  is  an  unknovm  but  non- 

random  object  field  in  aperture  R1  of  Figure  10,  EN(?2)  is 
a  zero-mean,  spatially  white  noise  field  representing  addi¬ 
tive  background  noise,  h21(r' ,r)  is  the  general  impulse  re¬ 
sponse,  and  hpS(r*,r)  is  the  free-space  impulse  response 
defined  by  Eq  (47).  Substituting  Eqs  (1J)  and  (2J)  into 
Eq  (59)  and  expanding  yields 


< IMSE  > 


■  /// //  ET(u)E*(m)hFe(v,u)h„q  (n,m) 

T>T“»T-»r-kr>^-  -**  ro  r  o 


R3R2R2R1R1 


x  hFS*(v,r3)hpS(n,r3)du  dm  dv  dn  d?3 


-  2  Re 


/////EI(5)Er(m)hFS(v';i)hFs'(','r3> 


R3R2R2R1R1 


x  <  h21*(n,m)h21(n,r3)>du 


u  din  dv  dn  dr,J 


-  2  Re 


II II ET(u)<E  *(n)>hFq(v,u] 

LR3R2R2Ri 


x  hFS#(v,r3)<h21(n,r3)>du  dv  dn  dr. 


+  1 1 1 1 1  ET(u)ET#(m)<h;>1  (v,u)h9  *(n,m) 

n  n  n  n  n  ^ 


R3R2R2R1R1 
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r 


*  —  —  —  —  —  —  — 

X  h21  (v,r3)h21(n,r3)>du  dm  dv  dn  dr. 


)><h21(v,u) 


■* 

-  2Re  J/JJe.IuXE  *<n)Xh <v,u 

lR3R2R2Rl 

« 

♦  -  —  _  _ 

x  h21  (v,r3)h21(n,r3)>du  dm  dv  dn  drg 
+  /  //<EN(^)EN*(H)><h2i*(^’?3) 


*  -  —  M  _ 

x  h21  (v,r3)h21(n,r3)>du 


R3R2R2 


x  h21(n,r3)>dv  dn  d?3  ;  u,  m€R1  v,  n< 


The  third  and  fifth  terms  of  Eq  (3J)  are  zero  because 
the  noise  field  is  zero-mean.  Since  the  additive  background 
noise  is  spatially  white  noise 


<EN(v)EN*(n)>  .  NoJ(v~H) 


For  propagation  through  the  turbulent  atmosphere,  the 
general  impulse  response,  h21(r',r),  of  Eq  (3J)  is  replaced 
by  the  turbulent  impulse  response,  hTB(r‘,r),  The  turbulent 
impulse  response  is  written  as 


hTB(r',?)  =  hFS(r' ,r)exp  [X(r' ,r)  +  j  $(?»,?)] 


Substituting  Eqs  (50)  and  (4J)  into  Eq  (3J)  yields 
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<IMSE>  =  f  f  f  f  f  ET*(m)hv~(v,u)hfrc*(n,m) 

o  n  n  n  n  A  r  ^ 


R3R2R2R1R1 


X  hFS*^*^3^hFS^»^3^<:‘|1  “  exp[X(n,m)] 
x  exp  [X(H,?3)]  exp[-j  0(n,m)  +  i(p(n,r3)] 

+  exp  [X(v,u)  +  X(v,r3)]  exp[j  $(v,u)] 
x  exp[- j  0(v,?3)]  +  exp  [X(v,u)  +  X(n,m)J 
x  exp  [X(n,?3)  +  X(v,?3)]  exp  [j  0(v  ,G)] 
x  exp[-j  0(v,r3)]exp[-j  0(nfm)  -j0(n,r3)] 


x  du  dm  dv  dn  d?3  + 


[no/(Az)2] 


//<  exp[2  X(n)]>dr 


R3R2 


2  <^r3 


;  u,  m€R1  v,  n€R2 


where  the  second  term  of  Eq  (3J)  is  expanded  using  the  iden- 

* 

tity,  Re(w)  =  (w+w  )/2  ,  and  a  change  of  variables 

_  _  _  _  # 

(u”*-*-m  and  v*«-*-n)  is  performed  on  the  w  term  of  the  expansion. 
Equation  (5J),  specifically  the  second  term,  can  be  simplified 
even  more.  Using  Eqs  (61)  and  (91),  it  can  be  shown  that 


( 
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<  exp  £2  X(n)]  >  =  1 


( 6J ) 


Substituting  Eq  (6J)  into  Eq  (5J)  yields 


<  IMSE  > 


-  /////  ET(u)E*(m)hFq(v,u)h  *(n,m) 

n  n  n  n  n  A  CO  CO 


R3R2R2R1R1 


x  hFS#(v,r3)hFS(n,r3)<|l  -  exp  [X(n,m)j 
x  exp  exp[-j  0(n,m)  +  j  0(n,r3)] 

+  exp  [X(v  ,u)  +  X(vfr3)]  exp  [j  0(v,u)] 
x  exp  [-j  0(v,?3)]  +  exp  [X(v»u)  +X(n,in)] 

x  exp  [X(n,r3)  +  X(v,r3)]  exp  0(v,u)J 

* 

x  exp  [-J  0(v,r3)]  exp  [-J  0<R  ,m)  +  j  0(H,?3)] 

< 

[no/(\z)2] 


> 


x  du  dm  dv  dn  dr3  + 


IU 


x  j  j  ar2  d?3  ;  u.  m€R1  v,  n€R2  (7J) 

R3R2 


Still  more  simplification  is  possible.  Since  Rg  is  a 
circular  aperture  of  radius  dg/2  and  R3  is  a  circular  aper¬ 
ture  of  radius  dj^/2  (same  radius  as  R^) ,  Eq  (7J)  reduces  to 
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<  IHSE  >  =  ///// Et(u)Et*(H 

n  n  r>  r»  n 


R3R2R2R1R1 


—  —  —  *  mm 

rn)hFS(v,u)hFS  (n,m) 


x  hFS  (v,r3)hpS(n,r3) <  |l  -  exp  [X(n,m)j 
x  exp  [X’(n,r3)]  exp[-j  0(n,m)  +  j  ^)(n,r3)] 

+  exp  [X(v  .u)  +  X(v,r3)]  exp  [i*  (v,u)] 
x  exp  [-j  0(v,?3)]  +  exp  [*<v  ,u)  +  X(n,m)] 

x  exp  [X(n,?3)  +  X(v,?3)j  exp  [j  0(v,u)j 

« 

x  exp[-j  0(v,r3)]  exp[-j  (f)(n,m)  +  j0(n,r3)]^> 


x  du  dm  dv  dn  dr~  +  N  D„ 

3  o  fo 


;  u,  m€R.  v,  n£R,  (60) 


where 


Df0  -[77V2/4H: 


(29) 


Equation  (60)  is  the  average  IMSE  of  the  image  field 
generated  by  the  CMF  Receiver  when  propagation  is  through 
the  turbulent  atmosphere  and  the  additive  background  noise 
is  zero-mean,  spatially  white  noise.  If  the  object  field 
is  also  isoplantic,  then  the  first  term  of  Eq  (60)  becomes 
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spatially  white  noise.  The  same  result  applies  however  for 
the  conditions  of  free-space  propagation  and  zero-mean,  spa¬ 
tially  white  noise.  This  is  shown  by  letting  X(r',r)  and 
0(r',r)  in  Eq  (5J)  equal  zero.  The  first  term  of  Eq  (5J) 
becomes  zero  and  Eq  (5J)  reduces  to 


<  IMSE  >  =  Nq/(Xz)2 


(10J) 


Since  R2  is  a  circular  aperture  of  radius  d2/2  and  R3 
is  a  circular  aperture  of  radius  d.^/2  (same  radius  as  R^ , 
Eq  (10J)  is  equivalent  to 


<  IMSE  >  =  N  D 

o  fo 


(61) 


where 


Df0  -  [nv2/4Xz]2 


(29) 
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Appendix  K  -  Sufficient  Condition  for  Minimizing  the  Average 
IMSE  of  the  Image  Field  of  the  Multiplicative  Phase  Receiver 

This  appendix  derives  Eq  (64).  Eq  (64)  is  the  sufficient 
condition  that  the  optimum  multiplicative  phase  of  the  MP 
Receiver  must  satisfy  for  the  average  IMSE  of  the  receiver 
image  field  t.o  be  minimized.  The  average  IMSE  of  the  image 
field  generated  bv  the  MP  Receiver  is 


< IMSE  > 


=  f<  [s(r3)  -  EI4p(r3)]  >dr. 


(63) 


where  E(r3)  is  the  field  for  free-space  propagation  and  is 
expressed  as 


E( 


r3^  =  J /EI(rl)hFS(r2,rl)h 
R2R1 

+  lE  N(? 


FS  ^2*^3^  d^l  d**2 


2)hFS  ^2,r3)d?2 


(IK) 


Using  Eqs  (23),  (51)  and  (53),  the  image  field  EMp(r3) 
generated  by  the  MP  Receiver  is  expressed  as 


EMP(*V 


J  J  EI(r1)exp[j  0  (r2)]  h21(?2,r1) 
(r2,r3)dr1  d?2  +  /  EN(?2) 


R2R1 


x  exp[j  0(?2)]hFS*(r2,?3)dr. 
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(2K) 


i 


( 


where  the  optimum  multiplicative  phase  0Q(r2)  has  been  re¬ 
placed  by  an  optimum  phase  estimate  0(?2).  In  Eqs  (IK)  and 
(2K),  Ej(r^)  311  unl<noWT1  hut  nonrandom  object  field  in 
aperture  R1  of  Figure  10,  EN(?2)  is  a  zero-mean,  spatially 
v/hite  noise  field  representing  additive  background  noise, 
h2i(r',r)  is  the  general  impulse  response,  and  hpS(r',r)  is 
the  free-spac-e  impulse  response  defined  by  Eq  (47).  Substi¬ 
tuting  Eqs  (IK)  and  (2K)  into  Eq  (63)  and  expanding  yields 


< IMSE > 


-  /////  EI(u)EI  (m 

T~*  T-l  T-»  "** 


R3R2R2R1R1 


_  —  *  _  _ 

)hFS(v,u)hps  (n,m) 


x  hpS  ( v,r3)hpS(n,r3)du  dm  dv  dn  d?3 


-  2Re 


)hpS(v,u)hpS*(v,r3) 


J  J  J  J  J*  i(u)EI*(m 

_  RgRgRgRj^Rj^ 

x  hpS(n,r3)<exp[-j  0(n)]  h21*(n,m)  > 
x  du  dm  dv  dn  d?3J 

J  J  J  J  EI(u)<EN*(n)><exp[-j  0(n)]> 

L R3R2R2R1 


-  2Re 


^  M  ^  mm  ^  mm 

x  hps(v,u)hps  (v,r3)hpg(n,r3)du  dv  dn  drg 


*  IJIJU 


R3R2R2R1R1 


j\ u)E j  (m)hpg  ( v i Tg ) hpg ( n | Tg ) 


x<exp[j0(v)  -  J  0(n)]h21(v,u)h21*(n,i) 
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I 


x  du  dm  dv  dn  dr. 


-  2  Re 


////.,« 

LR3R2R2R1 


)<En  (n)>hFS*(v,?3) 


x  hFS(n,r3)<exp[j  (f)(v)  -  j  0(n)]h21(v,u)> 


x  du  dv  dn  d?3 
♦  ///<  EN(v)EN*(H)>h 


R3R2R2 


iFs*(v,?3)hFs(H,?3) 


x  <  exp  [4  (v)  -  j  $m]  >dv  dn  dr. 


;  u,  m€R1  v,  n€R2  (3K) 


The  third  and  fifth  terms  of  Eq  (3K)  are  zero  because 

the  noise  field  is  zero-mean.  Now  let  the  optimum  phase 
A 


estimate,  (p  ( r2)  be  written  as 


f(?2)  -  0O(?2)  *  €0<?2> 


where  0Q(r2)  is  the  optimum  multiplicative  phase  and  €<^(r2> 
is  the  phase  error  of  the  optimum  phase  estimate.  Substitut¬ 
ing  Eq  (4K)  into  Eq  (3K)  yields 


< IMSE  > 


=  J  J  J  J  J  ET(U)E*(m)h 

R^R^R.R, R,  1 


FS(v,u)hps  (n,m) 


3  2  2  1  1 


r 


where  the  second  term  of  Eq  (3K)  is  expanded  using  the  iden¬ 
tity  Re(w)  =  (w+w  )/2  .  To  obtain  a  condition  for  0o(r2) 

such  that  the  average  IMSE  is  minimized,  take  the  derivative 
of  Eq  (5K)  with  respect  to  £,  let  €  =  0  ,  and  set  the 

result  equal  to  zero: 

i<IMSE>/^£|  =/j0„(n)  //// E  (u)E.*(m) 

'6-0  R2  «  R3RER1R1 

—  —  ♦  —  —  — 

x  hFS(v,u)hFS  (v,?3)hFS(n,r3) 

x  <exp[-j  0o(^)]h21  din  dv  d?3  dn 

mi 

+  J j0;(n)  -  J  J  J  J Et  (u)ET(m)h„„*(v,u) 

R2  L  R3R2R1R1 

x  hFS(v,r3)hFS*(n,r3)<exp  [j  0Q(n)j  h21(n,in) > 

x  du  din  dv  d?3  dn 

+  f  J&i  n)  J  J  f  J  Et(u)Et  (ffi)hF„(v,r3) 

R2  L  r3R2R1R1 

x  hFS*(n,r3)<exp[-j  0Q(v)  -  j  0Q(n)j  h21(n,u) 


(m)hFS(v,r3) 


x  h21*(v,m)>du  din  dv  d?3  J  dn 

+  J  J06(n)  -  J  J  J  J  Ei(0)Ei 

R2  L  R3R2R1R1 


)EI*(m)hFS*(v,r3) 


i 


hFs(n,r3)<exp  [j  -  j  0Q(n)  j  h21  ( v,u) 


*  —  —  _  _ 

x  h21  (n,m) >du  dm  dv  dr. 


+  /j  &(n)  J  J 


dn 


*  — 


R,  £  I  SX  ~~f*  - FS'->‘3' 

d  w  o  2 


<EN  (v)EH(n)>hFS(vfr3) 
x  hpS  (n,  r^X  exp  [-j  0q(  v)  +  j0o(n)]>dv  dr. 


+  J  3  0Jn)  -  J  J< 

R-  L  r3r2 


dn 


EN(v)EN*(H)>hFS*(v,?3) 


x  hpS(n,r3)<exp  [j  “  J0o(n)]>dv  dr; 


dn  =  0 


/ 

( 


;  u,  mCRp  v,  n€R2  (6K) 


where  the  fourth  term  of  Eq  (5K)  is  expanded  to  give  the 
third  and  fourth  terms  of  Eq  (6K),  and  the  fifth  term  of 
Eq  (5K)  is  expanded  to  give  the  fifth  and  sixth  terms  of 
Eq  (6K).  Equation  (6K)  is  true  if  the  sum  of  the  terms  in 
brackets  is  equal  to  zero.  Setting  the  sum  equal  to  zero 
and  performing  a  change  of  variables  (u-*-*-m)  in  the  third 
term  of  the  sum  yields 


IJJU 


R3R2R1R1 


*  —  —  * 

j'-JEj  (m)hFS(v,u)hFS  (v,?3)hFS(n,r3) 


x  <  exp  [- j  0o(n)]  h2i  (n,m)>du  dm  dv  dr. 
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r 


-  II IL 


rJJA*1  (u)EI(i)hFS*(^G)hFS(^'?3)hFS#(H*?3) 

o  2  1  1 


x  <  exp  [j  $Q(n)]  h21  (n,m)  >du  dm  dv  dr 
+  I f  J Jxx  (u)Ei(m)hFS(v,r3)hFS  (n,?3) 


^3^2 ^1^1 


*  <exp[-j0o(?)  +  J0o(F)]h21*(7,u)h21(H,5)> 


x  du  dm  dv  dr. 


~  1 1 1 1 E 

p  p  r  >  O  -L 


R3R2R1H1 


(S)hFS*(;,?3)hFs(H,?3) 


X  <exp[j0o(7)  -  J0o(H)]h21(7,C)h21*(B,K)> 


x  du  dm  dv  dr. 


+  //<EN  <v)EN(n)>hFS(v,r3)hFS*(n,r3) 

3  2 

x  <exp[-J0o(5)  +  J0o(n)]>dv 
*  / /<E „<v)E  *(E)>h„*(v,?,)h 


(v)En  (n)>hFg  (v,r3)hps(n,r3) 


x  <exp[j0o(v)  -  j  0o(H)]>dv  d?2  =  0 


;  u,  m  CR.,  v,  n  €R„ 


i 


Note  that  the  first  and  second  terms,  third  and  fourth 


terms,  and  fifth  and  sixth  terms  of  Eq  (7K)  are  conjugate 
pairs.  Therefore,  by  using  the  identity  (w-w*)  =  2ilm(w)  , 
Eq  (7K)  is  simplified  to 


J  J  J  f  E jCuJEj  (m)hFS*(v,r3)hFS(n,r3)<exp  [-J  <>0(n)] 

^  R3R2R1R1 

21*(n,m)  ^hFS(v,u)  -  exp  [j  $0(v)]h21(v,u)|  > 


x  du  din  dv  dr3  > 


-  Im 


//< 


*  - .  * - 


R3R2 


<En(v)En  (n)>hFS  (v,r3)hFS(n,r3) 
) 

x  <exp[j0Q(v)  -  j0o(n)]>dv  dr. 


>  =  0 


;  u,  mCRj^  v,  n  e  R2 


(8K) 


Since  the  additive  background  noise  is  spatially  white 

noise 


<V7)en 


*(n)>  =  N  d^v-n) 


(50) 


Substituting  Eq  (50)  into  Eq  (8K)  gives 


Im 


» 

'  ////E i(5>Ei*<!5>hFs‘(7'?3)hFS(;;-?; 


3  2  11 
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t 


» 

X  <exp[-j  (J)o(n)]  h21*(n,m)  hFS(v,u)  -  exp[j0Q(v)J 


The  second  term  on  the  left  side  of  Eq  (9K)  is  zero 
because  the  quantity  in  <  }  is  real.  Therefore,  Eq  (9K) 

is  simplified  to 


For  propagation  through  the  turbulent  atmosphere,  the 
general  impulse  response,  h21 (?•,?),  of  Eq  (10K)  is  replaced 
by  the  turbulent  impulse  response,  hTB(r',r).  The  turbulent 
impulse  response  is  written  as 

hTB(r',r)  =  hFS(r' ,r)exp  [X(r' ,r)  +  j0(r',r)]  (UK) 

Substituting  Eq  (11K)  into  Eq  (10K)  yields  Eq  (64): 


T 


Im 


//// 

R3R2R1R1 


—  *  _  —  * 

EjCuJEj  (m)hFS(v,u)hFS  (n,m)hFS  (v,?3) 


x  hFS(n,r3)<exp  [X(n,fn)J  exp[-j  (f)(n,m)  -  j$0(n)] 
r  jl  -  exp  [X(v,u)j  exp  [j  (j)  (v,u)  +  J$0(v)]j> 


x  du  dm  dv  dr_  >  =  0  ;  u,  m  €  R.  v,  n£R, 
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(64) 


Appendix  L  -  Minimum  Average  IMSE  of  the  Image  Field 


of  the  Multiplicative  Phase  Receiver 


This  appendix  derives  Eq  (65),  which  is  the  expression 
for  the  minimum  average  IMSE  of  the  image  field  generated  by 
the  MP  Receiver.  The  average  IMSE  of  the  image  field  gener¬ 
ated  by  the  MP  Receiver  is 


< IMSE  > 


=  J  J  f  f  J  Er(u)E*  (m)h„Av,u)h*  (n,m) 

R2  LR3R2R1Ri 

X  hFS*(v,r3)hFS(n,r3)du  dm  dv  d?3  j  dn 

-  I  Jiff  ET(u)ET  (m)hFq(vfu)hFq  (v,r„) 

rj  P  P  P  P  ’*■  *  ^  ^  ^ 

R2  LR3R2R1R1 

x  hFS(n,r3)<exp[-j  $0(n)  ]  h21*(n,in)  > 
x  du  dm  dv  d?3  j  dn 

-  /  III!  E*(u)ET(m)hFq  (v,u)hp„(v,?o) 

RP  P  R  P  ***  ^ 

X  hFS#(n,r3)<exp  [j  $Q(n)]  h21(n,m)> 
u  drn  dv  d?3  J  dn 

+  R  [/ RR  /EI(u)EI  ^m)hFS  (v»r3^hFS(n' 


x  du 


V 
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x  <exp[j  0o(v)  -  j  0o(n)]  h21(v,u)h21*(n,m)> 
x  du  dm  dv  d?3  J dn  +  N0Df0 

;  u,  m  €  1^  v,  n  €  R2 

(1L) 

Equation  (1L)  is  obtained  from  Eq  (3K)  in  the  following 

way.  The  third  and  fifth  terms  of  Eq  (3K)  are  zero  because 

the  noise  field  is  zero-mean.  The  second  term  of  Eq  (3K)  is 

* 

expanded  using  the  identity  Re(w)  =  (w+w  )/2  .  The  last 

term  of  Eq  (3K)  becomes  NQDfobecause  the  noise  field  is  spa¬ 
tially  white  noise.  This  is  shown  in  Appendix  J.  The  opti¬ 
mum  multiplicative  phase  estimate  0(?2)  is  replaced  by  the 
optimum  multiplicative  phase  0Q(r2) .  Finally,  the  integral 
order  is  rearranged  to  facilitate  a  forthcoming  substitution. 

The  optimum  multiplicative  phase  which  minimizes  the 
average  IMSE  given  by  Eq  (1L)  must  satisfy  the  sufficient 
condition  derived  in  Appendix  K: 
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In  its  present  form,  Eq  (10K)  cannot  be  easily  substi¬ 
tuted  into  Eq  (1L).  However,  using  the  identity 
Ira(w)  =  -i(w-w*)/2  ,  Eq  (10K)  is  put  in  the  more  useful 

form 


/// f  EI(u)E1*(m)hFS*(v,r3)hFS(n,r3) 

1 

[j0  (v)  -  j  0o(n)]h21(v,u)h21*(S,ni)>du  dm  dv  dr. 


R3R2R1R1 


x  <  exp 


=  /// / EI(u)EI*(m)hFS(v,u)hFS*(v,r3)hFS(n,r3) 


R3R2R1R1 


x  <exp[-j  0o(n)]h21*(n,S)>dn  dm  dv  dr3 

-  J  J  J  J  ^ t  (u)Ej(m)h-,g  (v,u)hpg(v,r3)hpg  (n,r^ 

R3R2R1R1 

x  <exp  [j4>0<s>Ki  (n,m)>du  dm  dv  d?3 

+  J  J  f  I  E*(u)EI(m)hFS(v,r3)hFS*(n,r3) 

R3R2R1R1 

x  <exp[-j0Q(v)  +  j  tf)0(n)]h21*(v,u)h21(n,m)> 


x  du  dm  dv  dr3  ;  u,  m  €Rj  v,  n€R2  (2L) 

Substituting  Eq  (2L)  into  the  fourth  term  of  Eq  (1L) 
yields 
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<IMSE>  =  /////  ET(u)ET*(m 

ri  T"»  r*»  n  ^ 


R3R2R2R1R1 


)hFs#(v,?3)hFS(n,?3) 


x  £hFg(v,u)hFS*(n,m)  -  2  hFS*(n,m)<  exp  [j  (f)o(  v)] 
x  h21(v,u)>  +  <exp[j  0Q(v)  -  j  <f)o(n)]h21(v,u) 
x  h21  (n,m)>j  du  dm  dv  dn  d?3  +  NoDfo 


;  u,  m€R1  v,  n€R0  (3L) 


where  a  change  of  variables  (u-*-*-m  and  was  performed 

on  the  third  term  of  Eq  (1L)  and  the  third  term  on  the  right 
side  of  Eq  (2L).  For  propagation  through  the  turbulent  atmos¬ 
phere,  the  general  impulse  response,  h21(r',r),  of  Eq  (3L)  is 
replaced  by  the  turbulent  impulse  response,  hTB(r',r).  The 
turbulent  impulse  response  is  written  as 


hTfi(r',r)  =  hFS(r' ,r)exp  [X(r*  ,r)  +  j$(r',r)| 


Substituting  Eq  (4L)  into  Eq  (3L)  yields  Eq  (65): 
<IMSE>  =  J  J I J I  EI(u)EI*(m)hFS(v,u)hFS*(n,m) 


R3R2R2R1R1 


x  exp  [j$><  v,u)  +  J0o(O)]>  +  <  exp  [  X(v,u)] 

x  exp  [X(n,ra)j  exp  [j  (/)  ( v,u)  +  j(f)Q(v)  -  j$(n,m)] 
x  exp j  0o(n)]  >  j  du  dm  dv  dn  d?3  +  N0Df0 

;  u,  v,  n€R2  (65) 
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Appe.idix  M  -  Average  Coherent  Transfer  Function 


of  the  Multiplicative  Phase  Receiver 


This  appendix  derives  Eq  (66),  which  is  the  expression 
for  the  average  CTF  of  the  MP  Receiver  when  it  is  operating 
so  as  to  cancel  turbulence-induced  phase  fluctuations.  The 
image  field  generated  by  the  MP  Receiver  is 


EMP^r3^ 


=  J E(r2)exp  [j  0  (r  )] 


0(r'2)JhFS  (r2,r3)dr2 


where  hFs^r2,r3^  is  the  free-space  impulse  response  defined 


by  Eq  (47): 


hFS(r2,r3)  =  exp [jk(z  +  J ?2-r3 | 2/2z]/ j X z 


The  field  E(r2)  of  Eq  (53)  is  the  input  to  the  MP  Re¬ 
ceiver.  By  Eq  (51)  it  is 


E(r2)  =  Eq(p2)  +  EN(r2) 


where  EQ(r2)  is  the  output  field  of  the  propagation  channel 
and  EN(r2)  is  a  zero-mean,  spatially  white  noise  field.  For 
propagation  through  the  turbulent  atmosphere,  the  field  EQ(r2) 
is  expressed  as 


,<?2>  *  /EI(?l)hTB(?2’?l)d?l 


7  -jn 


V  -Vf 


where  EI(?1)  is  an  unknown  but  nonrandom  object  field  in 
aperture  of  Figure  10  and  where  *1tb^2’*'i^  is  the  turbulent 
impulse  response  defined  by  Eq  (49): 

hTB^2*^l^  =  exp[jk(z  +  J  r2-r>1 1  2/2z)]/j  Xz 

x  exp  [X(r2,r1)  +  j  0(r2,?1)]  (49) 


Substituting  Eqs  (51)  and  (1M)  into  Eq  (53)  yields 


R2R1 


EMP^3^  =  ^  J^EI(?1)exp  [j  <^0(r2^]hTB^r2’^l^ 

R2R1 

x  hps  (?2,r3)dr1  d?2  +  J E^Cr^) 

R2 

x  exp  [j  0o(r2)]hFS  (r2,r3)d?2  (2M) 

Since  X(?',r)  and  (f)( r',r)  are  random  variables,  the 
turbulent  impulse  response  given  by  Eq  (49)  is  a  random  func¬ 
tion.  Therefore,  the  image  field  given  by  Eq  (2M)  is  also 
random.  On  the  average  however  the  field  is: 

<EMp(r3)>  =  //  EI(r1)<exp[j0o(?2)]hTB(r2,r1)> 


R2R1 


*  _  V  - 


x  hFS  (?2,r3)dri  dr2 


where  the  second  term  in  Eq  (3M)  is  dropped  because  the  noise 


.  4  ■,  y  ■  .j  • ;  fJi  ‘  _ 


field  is  zero-mean.  Substituting  Eqs  (47)  and  (49)  into 
Eq  (3M)  yields 


exp(  jkJ?3J  2/2z)<EjkIp(r3)  >  =  (1/Xz)2  </exp(  jk  j  J  2/2z) 

R1 

x  Ej(r^)Jexp  [-jk(r1-r3)  .?2/z]< exp  [<X'(?2  , ?1 )]  > 


x  <  exp  [j0(r2,r1)  -  j  (Jo(r2)]>dr2  d?x 


(4M) 


where  the  average  of  the  right  side  of  Eq  (4M)  has  been 
written  as  two  separate  averages  of  X(r',r)  and  0(r' ,r) 
since  these  are  independent  random  variables  (Ref  4:1374). 
Furthermore,  X(r' ,r)  and  (f)(  r',r)  are  Gaussian  random  varia¬ 
bles  (Ref  4:1374)  with  the  statistics  given  by  Eqs  (61),  (71) 
and  (81).  Using  Eqs  (61)  and  (91)  it  can  be  shown  that 


< exp  [X(r2,r1)]>  =  exp(-yr,2/2) 


( 5M) 


When  the  MP  Receiver  is  operating  so  as  to  cancel  turbu¬ 
lence-induced  phase  fluctuations,  then  0o(r2)  .  -0(r2)  , 

or  equivalently,  0O(?2)  =  “0(*;2,O)  *  Substituting  this 

choice  of  0O(^2)  and  Eq  (5M)  into  Eq  (4M)  yields 

exp(  jk|r3|  2/2z)<EJ(Ip(r3)  >  =  (l/Xz)2exp(-y^/2) 


X  /exp(  jk|rj  2/2z)E];(?;L)/exp[-jk(r1-r3)  .?2/z] 

R1  R2 

x  <  exp[j  0(r2,r1)  -  j  0  (?2 ,0)]  >d?2  dr^  (6M) 

Following  the  procedure  used  in  Appendix  I,  the  spatial 
average  on  the  right  side  of  Eq  (6M)  is  written  as 

<  exp  [j  0(r2,?1)  -  j0(r2,O)]>  =  exp[-D^(0,?1)/2]  ( 7M) 

where  D^(0,r1)  is  the  phase  structure  function  defined  by 
Eq  (181).  Substituting  Eq  (7M)  into  Eq  (6M)  yields 

exp(  jk|?3|  2/2z)<EIiIp(r3)>  =  (l/\z)2exp[-X?£ /2^ 

x  /exp(jkjr1|2/2z)EI(r1)exp  [-D^  (0,?1)/2] 

x  /  exp[-j2  TTCr^^-rg)  .?2/XzJdr2  drx  (8M) 

R2 

where  the  wave  number  k  in  the  inner  integral  is  replaced  by 
2  7T/X  .  It  is  shown  in  Appendix  I  that  the  inner  integral 
of  Eq  (8M)  can  be  written  as 

/exp[-J2  7T(r1-r3).r2/\z]dr2  =  dgX zJ^  TTj^-rgl dg/X z) 
R2 

x  (2/\r^r3\)  ( 9M) 


1  r*.»  1.-  ■  i 

-vi '■  :•  f 


Substituting  Eq  (9M)  into  Eq  (8M)  yields  Eq  (66): 


exp(  jk  |r3  |2/2z)<EMp(r3)  >  =  (d2/2  X  z)exp(-y^2/2) 

x  J  ex p(  jk  |r1 1  2/2z)EI(r1)exp[-D^(0,r1)/2j 
R1 

x  Jx(  7T|r1-r3|d2/\z)(l/  )d?1 
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(66) 


Appendix  N  -  Minimum  Average  IMSEs  of  Multiplicative 
Phase  Receiver  Image  Fields  for  Point 
Sources  and  Incoherent  Object  Fields 


This  appendix  shows  that  for  a  point  source  field  the 
sufficient  condition  for  minimizing  the  average  IMSE  of  the 
image  field  generated  by  the  MP  Receiver  is  satisfied  when 
the  optimum  multiplicative  phase  is  chosen  so  as  to  cancel 
the  turbulence-induced  phase  fluctuations.  This  statement 
is  qualified  by  saying  that  propagation  is  through  the  tur¬ 
bulent  atmosphere  and  the  additive  background  noise  is  zero- 
mean,  spatially  white  noise.  For  the  above  choice  of  the 
optimum  phase,  the  minimum  average  IMSE  of  the  image  field 
of  a  point  source  field  is  derived.  The  result  is  Eq  (70). 
The  process  of  showing  that  the  sufficient  condition  is  sat¬ 
isfied  and  then  calculating  the  minimum  average  IMSE  of  the 
image  field,  is  repeated  for  an  isoplanatic,  incoherent  ob¬ 
ject  field  with  a  real  and  even  object  intensity  distribution 
function.  The  minimum  average  IMSE  expression  for  this  case 
is  Eq  (71). 


When  propagation  is  through  the  turbulent  atmosphere 
and  the  additive  background  noise  is  zero-mean,  spatially 
white  noise,  the  sufficient  condition  that  the  optimum  multi¬ 
plicative  phase  must  satisfy  is 
t 


Im  < 


//// 

R3R2R1R1 


EjtujEj  (m)hFS(v,u)hFS  (n,m)hFS  (v,r3) 


v 
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1 


X  hpS(n,r3)<exp  [x(n,m)j  exp[-j  0  (n,m)  -  j  0Q(n)j 
x  -  exp  [X'(v.u)]  exp[j  0(vfu)  +  J0o(v)]j> 


x  du  dm  dv  dr3  r  =  0 


;  u,  m€R^  v,  n  €  R2  (64) 


Equation  (64)  is  derived  in  Appendix  K.  For  a  point 

source  field  with  a  complex  amplitude  Eq,  the  term 

_  *  _ 

Ej(u)Ej  (m)  in  Eq  (64)  is  written  as 


El(u)Ei*(m)  =  E0Eo*  cf( u)(f*(m) 


Substituting  Eq  (IN)  into  Eq  (64)  and  letting 
0  (r  )  =  -0(r?)  gives 


Im  i  (l/Xz)2EoEo*<//hFS*(v,r3)hFS(nfr3)r<exp  [X(n)]> 
{  R3R2  L 

-  <  exp  [X(n)  +  X(v)]>]  dv  d?3  1=0  ;  v,  n  €  Rg 


=  0  ;  v,  n€R, 


Using  Eqs  (61)  and  (91)  it  can  be  shown  that 


<exp[X(n)]>  =  exp 

<exp  [X(n)  +  X(v)]>  =  exp£-^-2J 


r 


>1 


Substituting  Eqs  (3N)  and  (4N)  into  Eq  (2N)  and  using 
Eq  (47)  to  expand  the  free-space  impulse  responses  yields 


(1/Xz)4E0Eq  £exp(-^/2)  -  exp(-^r  ;.]// 

]  23 

X  exp[-J2  TT(n  -v)  .?3/X  zjd?3  dv  >  =  0  ;  v,  n 


€.  R, 


(5N) 


where  the  order  of  integration  is  switched  and  the  wave 
number  k  in  the  impulse  response  is  replaced  by  2  7T/X  .  Note 
that  the  inner  integral  of  Eq  (5N)  is  a  two-dimensional  spa¬ 
tial  Fourier  transform  of  aperture  R3  where 

(n-v)/\z  =  ?  is  a  two-dimensional  spatial  frequency  vector. 
Since  R3  is  a  circular  aperture  of  radius  d^/2  (same  as  the 
radius  of  R^,  the  inner  integral  is  actually  a  Fourier-Bessel 
transform  of  circ(2 |r3|/d1)  where 

;  |?3|^V2 

;  otherwise  (6N) 


circ(2  r3J/d1)  = 


From  a  table  of  Fourier-Bessel  transforms,  the  transform 
of  Eq  (6N)  is 

[circ  ( 2 1  ?3 1  /d  i )  j  »  dj.XzJj.f  7T|n-v|d1/Xz)/2jn-vJ  (7n) 

Replacing  the  inner  integral  of  Eq  (5N)  with  Eq  (7N) 
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•  :'<•  /  • 


yields 


Im  «  (l/Xz)4E0E0*^exp(- v^2/2)  -  exp(-^)  j^Xz/2) 
x  f  (  7T|n-v|d1  /  Xz)(l/  Jn-v|  )dv  1=0  ;  v,  n 

R2  J 


2 

(8N) 


Since  J^Cx)  and  all  the  other  terms  on  the  left  side  of 
Eq  (9M)  are  real,  the  left  side  is  zero  and  the  sufficient 
condition  is  satisfied.  Therefore,  for  a  point  source  field 
and  the  conditions  stated  at  the  beginning  of  this  appendix, 
Eq  (64)  is  true  when  $0(r2)  “  -0(?2)  .  For  this  choice 

of  $0(r2)  and  the  same  conditions,  the  minimum  average  IMSE 
of  the  image  of  a  point  source  field  with  a  complex  amplitude 
Eq  is,  from  Eq  (65) 


< IMSE  > 


(1/Xz)  E0E0  JJ f  h?s  (v,r3)hpS(n,r3) 


R3R2R2 


£l  -  2  <  exp  [*<*)]>  +  <exp  [X(v)  +  X(n)]>j 


x  dv  dn  d?3  +  NoDfo  »  v,  n€R2  (9N) 


Substituting  Eqs  (3N)  and  (4N)  into  Eq  (9N)  and  using 
Eq  (47)  to  expand  the  free-space  impulse  responses  yields 

<IMSE>=  (1/Xz)4E0E0*  £l  -  2  exp(-^2/2)  +  exp ( - ) 
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(ION) 


J  J exp  [- J2  TC(.-r?/\  z)  .v  1  dv  / 
R3R2  J  R. 


x  exp  £-j2  7T(r3/X  z)  St  Jdn  d?3  +  NQD 


fo 


;  v,  n€R, 


where  the  wave  number  k  in  the  impulse  response  is  replaced 
by  2  7T7X.  Note  that  the  second  integral  is  a  two-dimensional 
spatial  Fourier  transform  of  circular  aperture  R2  where 
-r3/X  z  =  -?  is  a  two-dimensional  spatial  frequency  vector. 
Note  that  the  third  integral  is  a  two-dimensional  spatial 
Fourier  transform  of  the  same  aperture  v/here  r^/Xz  =  f  is 
a  two-dimensional  spatial  frequency  vector.  Therefore,  with 
the  proper  argument  changes  Eq  (7N)  can  be  used  twice  in 
Eq  (ION)  to  give  Eq  (70): 

<  IMSE  >  =  (d2/2Xz)2EoEQ*  1  -  2  exp(-^/2)  +  exp(-y£) 

X  /[Jl(Trd2|?3l/^z)]  (l/l?3l)2d?3 

R3 

+  NcDfo  (70) 


Before  the  above  process  is  repeated  for  an  incoherent 
object  field,  the  object  field  is  restricted  to  be  isoplanatic. 
For  this  additional  restriction,  Eq  (64)  becomes 


( 

i 

i 


! 

i 
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J  J  J  J  EI(u)EI*(fn)hFS(v,u)hFS#(n,m 


R3R2R1R1 


_  *  —  — 

m)hFS  (v,r3) 


X  hpS(n,r3)<exp  [x(n)]exp[-j  0(n)  -  j$Q(n)] 
x  ^1  -  exp  [X(v)]  exp  [j  0(v)  +  j  0o(v)]]> 


x  du  dm  dv  dr.,  >  =  0  ;  u,  m€.R,  v,  n  €  R, 


(UN) 


For  an  incoherent  object  field  with  real  and  even  object 

—  —  *  _ 

intensity  distribution  function  I(u),  the  term  E-^ujE^  (m) 

in  Eq  (12M)  is  written  as  a  spatial  average  over  the  object 
field 


<EI(u)EI*(in)>  =  I(u)  </(u-m) 


(12N) 


Substituting  Eq  (12N)  into  Eq  (11N)  and  letting 


0o(r2)  =  -^)(r2)  gives 


Im '  II  J"l(r1)hFS(v,r1)hFS*(n,r1)hFS#(v,r3)hFS(n,r3) 

R3R2R1 

x£<exp  [X(n)]>  -  <exp[X(n)  +  X(v)J>jdr1  dv  d?3  » 

4 


;  v,  n  €r. 


( 13N) 


Substituting  Eqs  (3N)  and  (4N)  into  Eq  (13N)  and  using 


Eq  (47)  to  expand  the  free-space  impulse  responses  yields 


Im  t  (1 


/Az)4[< 


exp(-^f/2)  -  exp(-^2) 


w. 


x  exp  [- j2  7T(n-v)  .?3/X  z]d?3  </l(r1)exp[-j2  7T*( v-n)  .i^/Az] 

R« 


x  dr^  dv 


>  -  0 


;  v,  n  €  R, 


( 14N) 


where  the  order  of  integration  is  switched  and  the  wave  number 
k  in  the  impulse  response  is  replaced  by  2  TT/A  .  Note  that 
the  second  integral  is  a  two-dimensional  spatial  Fourier 
transform  of  circular  aperture  R  where  (v-n)/Az  =  7  is  a 
two-dimensional  spatial  frequency  vector.  Therefore,  with 
the  proper  argument  change  Eq  (8M)  is  used  in  Eq  (14N)  to 
give 


Im  <  (1 


/Az)4[< 


exp(-Vl/2)  -  exp(-V5,f) 


j  (d2Az/2) 

J  Jr(.  7T|n-vjd  /Az)(l/|n-v|)J'l(r1 ) 

R2  R1 

x  exp  [-J2  7T(v-n)  .r^X  zjdr^  dv  |  =  0  ;  v,  n€! 


(15N) 


Note  that  the  inner  integral  of  Eq  (15N)  is  a  two- 
dimensional  spatial  Fourier  transform  of  the  object  intensity 
distribution  function  within  circular  aperture  R^,  Since  the 
function  is  real  and  even,  its  transform  is  also  real.  Since 
J^Cx)  and  all  other  terms  on  the  left  side  of  Eq  (15N)  are 
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real,  the  left  side  is  zero  and  the  sufficient  condition  is 
satisfied.  Therefore,  for  an  isoplanatic,  incoherent  object 
field  with  a  real  and  even  object  intensity  distribution  func¬ 
tion,  and  the  conditions  stated  at  the  beginning  of  this  ap¬ 
pendix,  Eq  (64)  is  true  when  d>0(?2> .  -0<  ?2)  .  For 

this  choice  of  0Q(r2)  and  the  same  conditions,  the  minimum 
average  IMSE  of  the  image  of  an  isoplanatic,  incoherent  ob¬ 
ject  field  with  a  real  and  even  object  intensity  distribution 
function  I(?1)  is,  from  Eq  (65) 


< IMSE  > 


IJIU 


R3R2R2R1 


I(r1)hFS(v,r1)hFS*(n,r1)hFS*(v,r3) 


x  hFS(n,r3)  1  -  2<exp[X(v)]>  +  <exp  [X(v)j 
x  exp  [X(n)]>]  dr^  dv  dn  dr*3  +  N0Dfo 


;  v,  n€R2  (16N) 


Substituting  Eqs  (3N)  and  (4N)  into  Eq  (16N)  and  using 
Eq  (47)  to  expand  the  free-space  impulse  responses  yields 


< IMSE  > 


(1/Az)4  [i  -  2  exp(-v^/2)  +  exp(-y^)J 

x  /  f  Hr^ </exp[-j27T(r1-r3).v/X  z]dv 
R3R1  R2 

x  ^exp[-J27T(r3-r1).n/Az]dn  drj  d?3  +  NQD 


;  v,  n€Rg 


(17N) 
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where  the  order  of  integration  is  switched  and  the  wave  number 
k  in  the  impulse  response  is  replaced  by  2  7T/X  .  Note  that 
the  third  integral  is  a  two-dimensional  spatial  Fourier  trans¬ 
form  of  circular  aperture  R2  where  (r^-r^/Az  =  -?  is  a 
two-dimensional  spatial  frequency  vector.  Note  that  the 
fourth  integral  is  a  two-dimensional  spatial  Fourier  trans¬ 
form  of  the  same  aperture  where  (r^r.^  )/X  z  =  f  is  a  two- 
dimensional  spatial  frequency  vector.  Therefore,  with  the 
proper  argument  changes  Eq  (7N)  is  used  twice  in  Eq  (17N)  to 
give  Eq  ( 71 ) : 


< IMSE  >  =  (d2/2\z)2  £  1  -  2  exp(-V^2/2)  +  exp(-^)J 

x  //I(?i)[Ji(7rd2lV?il/^z)]  (1/IV?il): 


R3R1 


x  drx  dr3  +  N0Df0 
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for  incoherent  object  fields  the  useable  decrees  of  freedom  are 
shown  to  varv  with  the  signal-to-noise  ratio  of  the  imaging  sys¬ 
tem,  The  SDatial  domain  and  freauency  domain  are  shown  to  be 
related  when  certain  conditions  are  satisfied,  one  of  which  is 
space-invariance . 

The  operation  of  two  ideal,  adaptive  imaging  receivers  (the 
Channel-Matched  Filter  and  Multiplicative-Phase  Receiver)  is  de¬ 
scribed  and  their  imaging  performance  is  compared  using  turbulent 
coherent  transfer  functions  and  minimum  average  integrated  mean 
square  error  expressions.  Compared  to  the  Channel-Matched  Filter 
Receiver,  the  Multiplicative-Phase  Receiver  is  shown  to  generate, 
in  general,  image  fields  with  larger  minimum  average,  integrated 
mean  square  errors. 
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